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ABSTRACT 

Records,  from  a  three  year  period  of  contract  stope  mining 
at  the  Beaverlodge  mine  of  Eldorado  Mining  and  Refining,  are  used 
to  develop  and  test  a  method  of  contract  system  analysis. 

This  method  produces  from  the  contract  or  production  records, 
reliable,  independent  and  unbiased  estimates  of  the  labour  and 
materials  content  of  the  contract  system  elements.  It  appears  that 
the  technique  would  be  of  value  to  any  industry  where  each  contract 
includes  mixtures  of  work  units. 
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INTRODUCTION 


This  thesis  is  concerned  with  the  use  of  multiple  regression 
techniques  to  analyze  the  underground  mine  production  data  from  the 
Beaver lodge  Operation  of  Eldorado  Mining  and  Refining.  The  purpose 
of  this  analysis  is  to  determine  a  basis  for  evaluation  of  con¬ 
tract  payments . 

The  Beaverlodge  Operation  of  Eldorado  Mining  and  Refining 
Limited  is  located  in  northwestern  Saskatchewan  and  produces 
uranium  oxide  from  1400  tons  per  day  of  ore  mined  almost  exclusively 
by  underground  methods.  Comprehensive  descriptions  of  the  property, 
geology,  mining  methods  and  production  processes  have  been  provided 
by  the  mine  staff.  (1,2)  A  longitudinal  section  and  surface  plan  of 
the  mine  is  shown  in  Fig.  1. 

The  mining  method  which  has  evolved  to  accommodate  the  high 
grade  irregular  orebodies  of  the  mine  may  be  classified  as 
selective  cut  and  fill.  Briefly,  this  consists  of  breaking  and  re¬ 
moving  successive  longitudinal  slices  from  the  underside  of  the  ore- 
body,  the  void  left  by  each  removal  being  filled  with  hydraulically- 
placed  mill  tailings. 

The  irregular  and  high  grade  nature  of  the  ore  forces 
repeated  departures  from  systematic  stoping.  Exploratory  test  holes, 
subdrifting,  slashing  and  raising  are  common  interruptions  to  the 
standard  break,  scrape  and  fill  cycle.  In  other  flatter  and  irregular 
orebodies,  problems  of  ore  removal  and  filling  cause  retention  of 
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several  slices  of  broken  ore  before  filling.  Thirty  to  forty  stopes, 
distributed  throughout  the  mine,  are  required  to  maintain  production. 

The  usual  operating  crew  for  a  stope  is  comprised  of  two 
miners  on  each  shift  working  on  piece-work  rates.  Mine  supervision 
and  the  geological  staff  control  the  grade,  mining  methods  and 
material  supply  by  daily  inspection,  but  the  workmen  are  mainly 
responsible  for  work  methods  and  material  usage. 

Piece-work  payments  are  administered  under  what  is  essentially 
a  sub-contract  system.  At  the  end  of  each  month  the  work  done  in 
the  month  is  measured  and  converted  to  dollars  by  means  of  straight-* 
line  rate  factors.  The  sumsof  these  converted  dollars  are  credited 
to  the  workmen  and  a  debit  is  made  for  total  hours  worked  and  also 
for  some  of  the  material  used.  The  credit  balance,  if  any,  is  paid 
in  dollars  to  the  workmen  in  proportion  to  the  hours  worked.  Records 
of  each  month's  measurements  and  calculations  for  all  stopes  have  been 
retained  since  the  beginning  of  mining  operations. 

Over  forty  different  types  of  work  are  recognized  in  the  con¬ 
tract  system.  Of  these,  thirty  three  which  comprise  99 %  of  the 
records  have  been  retained  for  analysis.  Some  nine  more  recorded 
values  such  as  material  used  and  time  spent  have  also  been  retained. 

At  first  glance  it  would  seem  that  contract  records  make 
analysis  of  stoping  production  a  simple  matter.  The  number  of  tons 
produced  divided  into  the  total  contract  payment  plus  basic  wages 
gives  an  overall  dollar  per  ton  cost  figure  for  each  stope.  Assign¬ 
ment  of  this  overall  cost  to  individual  types  of  work  done  in  the 
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stope  is  more  involved  but  can  be  estimated  by  proportioning  the  total 
cost  among  measurements  of  the  production,  according  to  work  classifi¬ 
cations,  using  the  rate  structure.  In  this  manner,  one  could  obtain 
a  unit  cost  per  ton  for  each  type  of  work.  These  unit  cost  figures 
with  their  dependency  on  unit  contract  rates  can,  however,  be  mis¬ 
leading.  If  the  individual  contract  rates  do  not  accurately  reflect 
individual  unit  labour  and  materials  content,  then  one  has  an  invalid 
analysis.  A  type  of  analysis  is  required  which  will  provide  a  valid 
estimate  of  the  labour  and  materials  content  of  production. 

Administration  of  the  piece-work  system  is  intimately  related 
to  stope  productivity  and  costs.  Maximum  productivity  and  minimum 
cost  can  be  maintained  only  by  a  set  of  just  and  equitable  contract 
rates.  Rates  are  only  equitable  when  they  reflect,  as  closely  as 
possible,  the  labour  and  controllable  material  content  of  all  pro¬ 
duction.  Rates  are  just  only  when  they  have  been  estimated  by  a 
valid  and  controllable  method.  A  further  problem  that  is  a  constant 
concern  of  administrators  is  the  assessment  of  modifications  in  the 
mining  methods.  Proper  evaluation  of  contemplated  modifications  can 
only  be  made  when  accurate  estimates  of  the  labour  and  material 
contents  of  present  methods  are  available. 

The  normal  method  for  estimating  labour  and  materials  cost  for 
each  type  of  job  is  by  time  and  method  study  of  actual  work.  This 
method  has  been  employed  extensively  at  Eldorado.  The  results  from 

t 

studies  in  single- job  production  work ,,  such  as  development,  have  been 
gratifying.  Results  from  studies  of  stope  work  also  have  shown  some 
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promise.  However,  with  such  a  large  number  of  job  types  within  a  wide 
spectrum  of  conditions,  development  of  timely  and  accurate  estimates 
of  unit  production  rates  has  been  extremely  difficult.  If  one  could 
obtain  such  estimates  rapidly  and  accurately  from  contract  records, 
ore  production  (which  is  the  central  concern  of  the  mine)  would  be 
much  more  amenable  to  administrative  analysis  and  control.  The  major 
concern  of  this  thesis  will  be  the  resolution  of  this  problem.  An 
attempt  will  be  made  to  estimate  from  contract  records  reliable  and 
accurate  rate  factors  of  labour  and  material  requirements  for 
individual  stoping  jobs  which  can  be  used  in  the  administration  and 
control  of  contract  stoping. 


... 


PREVIOUS  INDUSTRIAL  USE  OF  REGRESSION  ANALYSIS 


The  application  of  analytical  statistics  to  engineering  and 
management  problems  is  of  recent  origin.  Techniques  which  agriculture 
and  biology  have  used  for  years  have  just  recently  been  applied  in  the 
fields  of  engineering  and  management.  The  lack  of  problems  or  data 
has  not  been  the  major  cause  for  this  disparity.  The  opposite  has 
generally  been  the  case.  Such  a  large  quantity  of  data  is  generated 
by  modern  business  practices  that  only  the  advent  of  the  digital 
computer  has  allowed  effective  analysis  of  this  data.  This  condition 
probably  explains  why  a  text  by  Tippet  in  the  early  forties,  and  later 
re-edited  by  Lyle  (3)  has  not  received  wider  usage. 

Closer  to  the  present  (1960-1964)^  several  articles  (4,5,6) 
have  been  published  on  the  industrial  use  of  regression.  Campbell  (4) 
compares  the  tool  of  regression  analysis  to  that  of  a  microscope  in 
medicine.  With  it,  facts  hidden  in  the  mass  of  modern  business  data 
can  be  resolved  sufficiently  for  action  to  be  taken.  W.  B.  Glassford 
(5,  6)  gives  a  number  of  concrete  examples  where  regression  analyses 
were  used  to  construct  incentive  rate  systems.  Regardless  of  the 
title,  the  article's  claims  for  regression  are  modest  and  Mr.  Glassford' 
methods  of  hand  computation  can  be  easily  applied  in  problems  of  small 
dimensions . 

Extensions  of  regression  and  related  multivariate  techniques 
have,  with  the  advent  of  the  digital  computer,  become  quite  numerous 
in  the  mining  industry  (7,  8,  9,  10,  11).  All  these  applications  in 
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the  industry  to  date  have  been  mainly  concerned  with  sampling,  material 
properties  and  reactive  processes.  In  most  mining  operations,  materials 
and  processes  are  the  most  tractable  and  predictable  factors.  The  in¬ 
tractable,  unpredictable  and  often  highest  cost  items  are  usually  due 
to  labour.  Further  extension,  therefore,  of  proven  statistical 
techniques  for  analysis  of  mining  labour  production  data  should  be  of 
value  to  the  industry. 


THEORY 


The  major  objective  of  a  production  contract  system  is  the 
evaluation  of  the  skill  and  effort  which  workmen  employ  in  the  pro¬ 
ductive  process.  Monetary  evaluation,  given  for  skill  and  effort,  is 
returned  to  the  workmen  as  incentive  pay  or  bonus.  The  usual  method 
of  evaluation  is  based  on  a  measurement  system  under  classified  con¬ 
ditions  for  a  set  period  of  time.  Units  of  production  measured  are 
converted  to  monetary  credit^ and  time  worked  is  debited^by  a  function 
operating  on  the  measured  values.  Workmen  are  paid  the  positive 
balances  while  negative  balances  are  absorbed  by  the  employer. 

The  method  of  conversion  of  units  of  work  to  units  of  value 

a 

used  at  the  Beaverlodge  Operation  of  Eldorado  Mining  and  Refining  up 
to  1964  was  a  simple  linear  equation.  Rates,  in  dollars  per  unit 
produced,  were  constant  throughout  the  range  of  production,  type  of 
mining  and  location  of  the  working  place  in  the  mine.  The  method  of 
debiting  time  worked  has  also  been  constant  and  linear.  (A  very  small 
amount  of  labour  has  been  charged  at  $1.40  per  hour  but  the  major  part 
has  been  charged  at  a  constant  rate  of  $1.56  per  hour.) 

This  system^ with  its  mine-wide,  time-constant  application, 
implies  that  productivity  is  constant  in  time  and  place.  The  linear 
nature  of  the  function  implies  that  production  units  are  independent 
both  of  each  other  and  of  the  level  of  production.  It  would  be 
desirable  to  know  if  these  implications  are  true^ and  if  the  "best" 
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functions  have  been  used  for  the  desired  evaluation.  To  solve  this 
problem  one  must  first  define  what  is  meant  by  "best". 

When  measuring  a  variable  (the  skill  and  effort  of  workmen 
certainly  is  variable)  by  use  of  one  or  more  indicators,  the  usual 
method  is  to  observe  simultaneous  occurrences  of  the  dependent  variable 
and  its  proposed  indicators  (independent  variables).  From  these 
observations  one  constructs  mathematically,  graphically,  or  intuitively 
a  function  which  operates  on  the  independent  variables  to  predict  (fit) 
the  dependent  variable.  The  differences  between  the  predicted  values 
and  the  observed  values  of  the  dependent  variable  within  the  set  of 
observations  are  now  used  to  evaluate  the  constructed  function.  One 
may  use  several  separate  criteria  of  evaluation  such  as  the  minimum 
sum  of  the  absolute  value  of  this  difference  or  residual,  or 
minimization  of  the  maximum  residual  value.  The  criterion  however 
which  is  the  most  tractable  mathematically  and  generally  used  is  that 
the  sums  of  the  squares  of  the  residuals  shall  be  a  minimum. 

Throughout  this  work  the  minimum  sum  of  squared  residuals;  i.e.  the 
least  squares  criterion,  will  be  considered  the  "best"  method  of  "fit". 

The  residual  value  from  a  least  squares  "fit"  of  function  to 
observations  has  a  second  and  very  important  statistical  use.  If  the 
frequency  distribution  of  the  residual  values  can  be  considered  normal 
with  a  mean  of  zero  and  independent  variance  then  one  has  exact  methods 
for  attaching  probability  statements  to  the  estimates. 
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The  contract  system  of  computing  payments  and  costs,  which  is 
a  linear  function  used  to  measure  or  predict  the  skill  and  effort  of 
workmen,  can  be  compared  to  any  other  function  applied  to  the  same  set 
of  observations.  The  comparison  can  be  made  at  two  levels,  the  form 
of  the  functions^ and,  if  the  functions  are  the  same,  the  parameters 
within  the  functions. 

Let  it  first  be  assumed  that  the  forms  of  the  functions  to  be 
compared  are  identical.  A  function  containing  the  observed  independent 
variables  (production  unit  measurements)  of  the  contract  system  will 
be  fitted  to  the  observed  hours  worked  (dependent  variable) .  If  the 
estimated  function  is  able  to  predict  the  dependent  variable  with  a 
significant  reduction  of  its  variance  (analysis  of  variance  tests) 
then  we  shall  conclude  that  the  form  of  the  function  is  valid.  One 
can  next  compare  the  parameters  of  the  function  with  those  given  as 
contract  rates.  If  the  absolute  difference  between  the  given  rate 
and  that  estimated  is  greater  than  that  which  could  be  attributed  to 
one  chance  in  twenty  (957o  confidence  limit)  then  one  can  say  that  the 

L 

given  contract  rates  are  biased  or  not  "best". 

The  above  assumes  that  the  residuals  will  not  be  significantly 
different  from  the  foregoing  normal  and  independent  specification. 

If  this  criterion  is  not  met  then  the  probability  statements  on  form 
and  parameters  will  lose  a  measure  of  exactitude,  the  amount  of  loss 
being  somewhat  determinable  by  the  level  of  the  statistics  used  to 


test  the  normality. 
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The  problem  can  be  formulated  using  two  differing  but  related 
models.  The  direct  approach  using  a  form  of  function  like  the  given 
system  may  be  termed  "Labour  Content  Estimation".  The  raw  estimates 
of  the  parameters^ in  units  of  hours  per  unit  of  production,  are  analagous 
to  the  type  of  estimates  made  from  time  and  methods  study,  the  usual 
method  for  constructing  a  contract  rating  system.  A  second  related 
method  or  model  attacks  the  problem  indirectly. 

The  indirect  method  may  be  termed  "Correction  Factor  Estimation". 
If  the  contract  system's  intent  is  to  measure  the  skill  and  effort  of 
workmen  in  dollars,  then  the  variations  in  this  dollar  total  for  equal 
periods  of  time  should  be  an  unbiased  minimum.  In  other  words  the 
given  criteria  of  the  "best"  could  directly  evaluate  the  rates  used  to 
calculate  incentive  earnings.  Let  us  consider  that  contract  rates, 
for  production  unit  values,  are  scaling  or  conversion  factors  to 
change  all  measurement  units  to  a  common  dimension  of  dollars.  Let 
the  mean  incentive  rate  plus  the  standard  contract  wage  be  the  contract 
time-conversion  factor.  Now  estimate  a  function  which  will  "fit"  the 
dollar  earnings,  as  represented  by  rate  conversion,  to  the  observed 
earnings  as  represented  by  observed  contract  hours  times  mean  earnings 
rate.  If  the  estimated  function  is  able  to  significantly  reduce  the 
variance  of  the  earnings  per  contract  an  "improved"  contract  rate 
function  has  been  constructed.  The  remaining  variance  in  earnings 
per  contract  may  be  assigned  to  skill  and  effort  of  the  workmen  and 
numerous  other  factors  which  might  affect  work  rates.  The  function 
has  only  removed  variance  in  earnings  due  to  inequalities  in  rates 
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and  in  this  sense  has  estimated  "improved"  rather  than  "best"  rates. 

If  the  estimated  values  of  the  coefficients  of  the  function 
are  significantly  different  from  1.00  then  we  can  apply  these 
parameters  as  multiplicative  factors  for  the  past  contract  rates  to 
produce  what  would  be  "improved"  rates.  Another  decided  advantage  of 
the  coefficients  estimated  from  this  model  is  that  all  estimates 
are  non-dimensional  and  can  be  compared  with  one  another. 

Both  of  the  foregoing  models  can  be  represented  mathematically 

by: 

Y  =  Bo  Xo  +  Bl  Xx  +-...  +  Bk  Xk  (1) 

In  the  Labour  Content  Estimation,  Model  I,  Y  is  the  total  hours 
worked  in  a  given  contract  period  in  a  given  working  place,  X^ ,  , . . . 

X  are  the  measured  values  of  the  "k"  different  types  of  production  in 
the  working  place,  (X  =  1)  and  B  ,  B...B,  are  the  rates  in  hours  per 

O  OK 

unit  of  the  "k"  different  units.  In  the  Correction  Factor  Estimation 
Model  II,  Y  is  the  contract  dollar  credit  in  a  given  contract  period 
in  a  given  working  place,  (XQ  =  1)  X-^,  X2...Xjc  are  the  calculated 
dollar  values  of  the  "k"  different  units  produced  in  the  working  place, 
and  Bq,  B^ ,  ^...B^  are  non-dimensional  factors  which  qualify  the  con¬ 
tract  rates  used  to  calculate  the  dollars  values  of  X^. 

The  calculation  procedure  to  estimate  the  parameters  of  a  linear 
model  fitted  to  observations  by  the  least  squares  criteria  can  be 
illustrated  for  a  three  dimensional  case.  Extension  for  the  "k" 


dimensions  of  the  model,  equation  (1),  can  be  made. 
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Let  the  model  be  Y=  BQ  XQ  +  B^  X^  +  B2  X2  (2) 

where  XQ  =  1  and  Xp  X2 ,  are  the  independent  variables  and  Y  the 

1 

dependent  variable.  One  has  N  observations  (N  ^  3)  and  wishes  to 
estimate,  from  this  set  of  observations,  coefficients  BQ,  B^  and  B2, 
which  will  give  the  equation  of  a  plane.  This  plane  is  to  be  at 
that  position  where  the  sums  of  squares  of  the  differences  between  the 
predicted  dependent  variable  and  the  observations  on  that  variable 
are  a  minimum.  The  usual  method,  partial  differentiation,  for  find¬ 
ing  the  minimum  of  a  function  is  applied  to  the  observations 
according  to  equation  (2)  and  results  in  the  following  simultaneous 
equations  called  normal  equations. 

S  (Y)  =  b0  N  +  bx  S  (X^  +  b2S  <X2)  (3) 

S  (YXp  =  bQ  S  (XL)  +  bx  S  (Xx2)  +  b2  S  X2)  (4) 

•. 

S  (YX2)  =  bo  S  (X2)  »  bj_  S  (Xj_  X2)  b2  S  (X22)  ,  (5) 

where  3  indicates  summation  over  N  observations.  Solving  by 


where  D  is  the  determinant  of  the  square  symmetric  matrix 
formed  by  the  sums,  sums  of  squares  and  sums  of  cross  products  of  the 
right  hand  side  of  equations  3,  4,  and  5^and  Dp  D2  and  are  similar 
determinants  formed  by  replacing  rows  1,  2  or  3  of  D  with  the  left  hand 
sides  of  equations  3,  4,  and  5. 

It  is  easily  seen  that  extension  of  the  above  procedure  to  "k" 
dimensions  makes  machine  computation  mandatory.  (For  this  study  there 
were  1275  sets  of  data  with  42  variates  per  observation) .  A  computer 
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program  prepared  by  Dr.  K.  Smillie  of  the  Department  of  Computing 
Science,  University  of  Alberta  (12)  computes  the  required  sums,  sums  of 
squares ) and  sums  of  cross  products  from  all  given  observations  and 
transformations  made  on  these  variables.  From  this  large  matrix  a 
subset  is  selected  as  independent  variables  for  regression  against 
a  chosen  dependent  variable.  The  regression  is  done  in  a  step-wise 
additive  manner  with  the  calculation  of  analysis  of  variance  between 
each  additional  term  added.  Control  of  the  order  of  the  regression 
calculation  can  be  external  or  may  be  left  to  the  computer.  Computer 
control,  which  selects  the  variables  in  the  order  in  which  they  will 
reduce  the  percent  variance  in  the  dependent  variable,  was  used 
throughout  the  following  analysis.  This  order  is  not  necessarily 
related  to  the  significance  of  individual  coefficient  estimates  but 
in  this  case  the  trends  are  generally  the  same. 

Tests  of  significance  on  individual  coefficient  estimates 
can  be  made  from  the  standard  error  of  estimate  calculated  for  the 
coefficient.  Other  statistics  which  can  be  used  to  evaluate  various 
conditions  of  the  regression  and  the  residual,  after  all  terms  have 
been  added,  are  produced  by  the  program  used  (12) . 

Their  use,  where  applicable,  is  demonstrated  in  following 


sections . 
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DATA  -  SOURCES,  PREPARATION  AND  PROCESSING 


In  most  experimental  work  the  apparatus  used  to  generate  the 
data  is  of  primary  importance  and  the  tools  of  analysis  receive  little 
attention.  For  this  study  of  mine  contracts,  controls  on  the  genera¬ 
tion  of  data  were  nonexistent  and  the  analysis  of  data  is  the  main 
concern.  For  this  reason^  then,  the  methods  used  to  sample  the  data  and 
prepare  it  for  analysis  are  examined  in  some  detail. 

The  sample  of  1275  contract  calculations  includes  over  357,  of 
all  stope  contract  records  for  the  period  Feb.  1962  to  Nov.  1964. 

The  contract  system  under  which  these  records  were  generated  is  given 
under  "Stoping  Contract  Regulations"  (Appendix  3,  pages  1  and  2). 

The  terminology  for  various  types  of  stope  work  is  fairly  common 
throughout  the  industry.  Detailed  descriptions  of  the  stoping  systems 
and  related  terms  are  available  in  a  comprehensive  treatment  of  the 
mining  operation  to  which  reference  has  been  made  (2) . 

The  numbering  system  under  which  the  various  contract  elements 
were  coded  for  computation  is  given  in  Table  1  (Appendix  3,  page 
4) .  It  will  be  noted  that  the  number  of  contract  elements  has 

been  reduced  somewhat  from  that  provided  by  the  contract  regulations. 

This  reduction  was  made  by  some  grouping  and  some  exclusion  during 
preparation  of  the  records  for  key  punching.  Notes  taken  during  this 
preparation,  given  in  Appendix  3,  page  5,  describe  the  nature  of  this 
grouping  and  exclusion.  Additional  information  for  classification  of 

the  data  was  obtained  from  the  mine  staff.  A  brief  description  of 


-2 


15 


classification  schemes  is  given  in  Table  2  (Appendix  3,  pages  6  and  7) 
The  enumeration  of  the  records  in  each  classification  follows 
(Appendix  3,  pages  8,  9  and  10). 

The  main  computer  program  used  in  the  stepwise  multiple  re¬ 
gression  was  a  library  program  of  the  Department  of  Computing  Science, 
written  for  Dr.  K.  Smillie,  by  Mrs.  W.B.  Payne  (12).  This  program  is 
a  general  stepwise  multiple  regression  program  written  in  FORTRAN  IV 
for  the  IBM  7040  which  has  32,000  words  of  core  storage.  Two 
magnetic  tape  units  are  required  for  intermediate  results,,  and  cal¬ 
culations  are  in  single  precision  arithmetic.  A  comprehensive 
description  of  the  program  is  available  in  the  Department  of  Computing 
Science,  University  of  Alberta,  Edmonton,  Alberta.  A  key  to  the 
portions  of  this  program's  output  retained  for  analysis  is  given  in 
Table  3  (Appendix  3,  pages  11,  12,  13  and  14). 

A  second  major  program  used  was  the  "Autoplot"  program  (13). 
This  program  used  the  IBM  7040  system  to  produce  a  typewriter  plot  of 
data.  The  7040  produces  punch  cards  which  are  manually  fed  to  a  card- 
reader/document-typewriter  system  to  produce  the  final  plot.  The 
program  is  written  in  the  Assembly  language  of  the  7040  system  and 
is  under  direct  control  of  the  system  monitor. 

Specific  programs  were  written  to  arrange  raw  data  and  the 
multiple  regression  output  in  a  form  adaptable  to  the  Autoplot  program 
As  these  small  programs  were  essentially  simple  data  handling  opera¬ 
tions  particular  to  the  immediate  problem  no  attempt  has  been  made  to 
record  and  maintain  them.  Descriptions  of  their  operation  given  in 
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following  sections  should  be  quite  sufficient  to  allow  one  to  write 
similar  FORTRAN  programs. 


MULTIPLE  REGRESSION  PROGRAMS  ON  DATA  AND 


DATA  CLASSIFICATIONS 


(A)  General  Procedure 

Experimental  procedures  for  data  analysis  follow  much  the  same 
path  as  those  in  any  other  course  of  experimentation.  A  sample  of  the 
available  data  is  tested  against  one  or  more  working  hypotheses.  The 
results  are  examined  and  experimentation  continued  with  larger  samples 
and  mo«*e  precise  hypotheses.  In  conclusion  a  general  function  of  the 
variables  is  proposed  and  tested  by  observing  the  reliability  and 
veracity  of  the  function  under  specific  conditions. 

The  first  tests  were  made  by  graphical  representation  of  the 
data.  These  tests  indicated  little  more  than  the  complexity  of  the 
problem  and  were  followed  immediately  by  regression  analyses  according 
to  the  mathematical  models  given  in  Theory. 

The  first  regression  tests  on  a  sample  of  about  50  per  cent  of 
the  observations  were  made  using  the  Labour  Content  Model.  The  re¬ 
sults  of  this  test  (14)  confirmed  the  hypotheses  that  multiple  linear 
regression  on  contract  records  would  produce  reasonable  measures  of 
the  labour  content  of  individual  production  units.  The  remainder  of 
the  data  were  gathered  and  subjected  to  a  similar  computation.  Re¬ 
sults  of  this  computation  are  given  in  Appendix  1  as  Program  3. 

Program  3  was  now  examined  in  detail  and  the  relationships  of 
the  estimated  labour  content  coefficients  and  the  rates  of  the  con¬ 
tract  system  established.  A  full  scale  test  was  now  made  using  the 
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Correction  Factor  Model.  The  output  of  this  test,  Program  60,  is 
given  in  Appendix  2.  Comparison  of  the  results  of  these  two  tests 
was  made  to  establish  the  relationship  between  the  two  methods  of 
approach.  As  was  expected,  a  direct  relationship  was  demonstrated. 

The  Correction  Factor  method  was  used  for  further  analysis  of 
the  data  according  to  several  different  classification  systems.  The 
classification  systems  were  as  follows: 

1.  Mining  Methods  -  6  classes  as  outlined  in 
Table  2  plus  a  seventh  class  made  by  divid¬ 
ing  B  type  stopes  into  Fay  and  Verna  areas, 
the  two  major  operating  divisions  of  the  mine. 

2.  Incentive  Pay  Levels  -  all  data  divided  into 
$0.50  per  hour  increments  of  incentive  pay 
rate . 

3.  Time  Series  with  a  6  Month  Moving  Average  - 
all  data  were  grouped  into  6  month  periods 
starting  February,  1962.  Three  sets  of  groups 
were  made  so  that  calculation  periods  were  two 
months  apart. 

The  results  from  all  the  above  classifications  were  compared 
with  the  results  from  a  single  regression  on  all  data  (Program  60) 
by  covariance  analysis.  This  is  a  method  in  which  one  compares  the 
overall  variance  in  the  predicted  value  using  a  set  of  classifica¬ 
tions  with  the  variance  of  prediction,  using  the  same  model,  on  the 
complete  data.  The  results  of  those  tests,  given  in  Tables  6,  7, 
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and  8  ,  Appendix  4,  indicated  quite  clearly  that  all  three  types  of 
classifications  were  valid. 

A  method  was  now  required  to  compare  the  parameter  estimated 
in  one  class  with  estimates  of  the  same  parameters  in  the  other 
classes  within  the  same  classification.  Exact  statistical  methods 
were  contemplated  for  this  task,  but  a  method  closer  to  industrial 
practices,  the  control  chart,  appeared  to  serve  the  purpose  well  and 
would  have  a  more  general  appeal.  FORTRAN  programs  were  written  to 
take  the  regression  results  of  a  set  of  classifications  and^using  the 
estimated  standard  deviation  of  individual  coefficient  estimates,  to 
construct  a  table  of  values  representing  the  95  per  cent  confidence 
limits  of  this  estimate.  (As  the  number  of  observations  in  any  one 
class  is  greater  than  70^ exact  normal  theory  was  used  to  give  95  per 
cent  limits  at  1„96  times  the  standard  error  of  estimate). 

The  parameter  estimates  plotted  have  also  been  altered  slightly 
by  a  weighting  function  to  make  them  more  comparable  to  the  contract 
rate  system  and  each  other. 

In  the  contract  system  the  rates  were  calculated  on  a  full 
8-hour  day  basis,  a  part  of  which,  due  to  travel,  lunch,  and  rest  periods 
was  non-productive.  The  regression  model  reflects  somewhat  this 

non . productive  time  by  inclusion  of  the  B0  parameter.  Now  since 

contract  rates  have  been  upward  adjusted  to  allow  for  the  unavoid¬ 
able  inclusion  of  non-productive  time,  it  seems  reasonable  that 
parameters  estimated  from  regression  should  be  treated  similarly. 
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A  multiplicative  factor  consisting  of  34  was  used  to 

*34  -  B0 

adjust  all  estimated  coefficients. 

(X34  represents  the  mean  value  of  the  dependent  variable  in  the 
regression  and  B0  the  coefficient  estimated  for  XQ  in  the  same 
regression) . 

After  examination  of  the  results  from  the  regression  pro¬ 
grams  for  evaluation  of  labour  variables  it  appeared  that  some 
useable  information  on  material  consumption  might  be  extracted 
from  the  contract  records.  As  prior  work  had  indicated  that  all 
classifications  were  valid  and  the  plotting  method  gave  a  concise 
illustration  for  evaluation,  it  was  decided  to  proceed  directly  to 
this  step  in  the  analysis.  Five  separate  sets  of  regression 
analyses  were  made  under  six  classifications.  The  computer  outputs 
are  on  file  (22)  and  the  plots  of  significant  parameters  are  given 
in  Figures  52-66. 

(B)  General  Tests  of  Methods 

(1)  Graphical  Representation  and  Analysis  of  the  Data 
For  much  experimental  work  the  plotting  of  data  in  graphi¬ 
cal  form  allows  one  to  infer  quite  useful  and  testable  functions  for 
further  analysis.  In  this  mining  contract  situation  the  dimensions 
of  the  data, namely  30  to  40  different  jobs  or  types  of  work,  make 
direct  graphical  work  impossible.  To  reduce  the  dimensions  somewhat 
a  number  of  the  variables  were  combined.  This  was  done  by  a  small 
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FORTRAN  program  reading  the  raw  data  from  magnetic  tape.  The 
combined  var iablesj-subdri f ting ,  ore  breaking,  slushing,  slusher 
floor,  millhole  construction  and  tramming  are  the  sums  of  the 
individual  units  in  each  class  as  grouped  in  any  one  observation 
(contract  period) .  (Table  1  gives  a  listing  of  these  variables 
under  their  respective  grouping  nomenclature) .  These  six  variables 
and  two  more,  the  number  of  blast  holes  drilled  and  the  pounds  of 
explosive  reported,  were  plotted  against  contract  hours  per  contract 
per  month.  Four  more  plots  were  also  made  of  development  footage 
and  tons  of  ore  broken  versus  blast  holes  drilled  and  pounds  of 
explosive  reported.  The  above  12  plots  are  given  in  Figures  2-13. 

Examination  of  these  plots  indicates  that  graphical  analysis 
is  of  little  quantitative  value  in  this  problem.  There  is  an  indi¬ 
cation  that  all  pairs  of  values  plotted  may  have  positive  linear  re¬ 
lationships  but  it  would  be  rather  unwise  to  draw  lines  to  represent 
those  relationships.  The  best  that  one  can  say  of  this  type  of 
analysis  is  that  it  clearly  indicates  the  requirements  of  more 

sophisticated  techniques  for  any  useful  solution  to  the  problem. 

(2)  General  Tests  on  the  Multiple  Regression  Method 

A  major  assumption  made  to  allow  the  use  of  regression  methods 
concerns  the  frequency  distribution  of  the  residual.  If  the  residual^ 
the  difference  between  the  observed  and  predicted  value  of  the 
dependent  variable*  has  a  frequency  distribution  which  is  random  and 
independent  than  one  may  validly  apply  statistical  tests  for  evalua¬ 
tion  of  the  regression  model  and  its  associated  parameters. 


. 

l  zb  i)oq  Oi'.i  b 

{ 

...  ..  • 

,n«  j  'tb 


22 


To  allow  one  to  visualize  the  relationship  of  predictions 
to  observations  and  of  observations  to  residualg,  ten  per  cent  of 
the  samples  of  the  observed  versus  predicted  values  and  residuals 
versus  the  observations  were  plotted  (Figures  14  and  15) ,  using 
the  Autoplot  system  (13) .  The  values  used  for  these  plots  are 
given  in  the  Reference  Data  Vol.  1  (15). 

Figure  14»  of  observed  versus  predicted  values,,  indicates 
that  while  the  prediction  equation  may  be  unbiased  (since  the 
number  of  points  above  and  below  the  equality  line  are  substan¬ 
tially  the  same)  the  grouping  tendency  around  the  value  of  300 
hours  per  month  leaves  some  doubt  as  to  randomness.  The  same 
grouping  tendency  is  also  evident  in  Figure  15,  observed  versus 
residual.  The  grouping  tendency  in  both  plots  may  be  ascribed 
to  the  fixed  time  period  for  contract  calculations  of  one  month 
and  the  relatively  constant  crew  strength  of  two  men.  True  non¬ 
randomness  of  the  residuals  can  be  assessed  by  consideration  of 
the  distribution  of  the  points  of  Figure  15.  Because  there  is  no 
apparent  tendency  for  the  values  of  the  residual  to  concentrate 
about  one  or  more  horizontal  lines  or  bands,  one  can  assert  that 
there  is  no  evidence  that  the  residual  is  not  random. 

Figures  14  and  15  show  no  strong  evidence  that  the  residual 
is  not  independent  of  the  dependent  variable.  Evidences  of  depend¬ 
ency  of  the  residual  are  usually  shown  in  these  types  of  plots  by 
fan  shaped  or  parabolic  trends  of  the  points. 
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For  more  exact  tests  on  the  residuals  the  summary  statistics 
of  these  residuals  which  complete  all  of  Dr.  Smillie's  programs 
may  be  considered.  These  statistics  (see  Appendix  2,  page  25) 
have  a  variety  of  uses, some  of  which  are  not  applicable  to  this  case. 
The  statistics  applicable  here  are: 

1.  Sums,  sums  of  squares,  and  mean  of  the  residuals  * 

2.  Coefficient  of  skewness  and  coefficient  of  kurtosis 
of  the  residuals  and  their  associated  standard 
deviations . 

The  remaining  statistics  may  be  used  to  test  the  serial 
nature  of  data.  These  might  be  used  to  indicate  that  the  proposed 
time  series  and  other  classifications  of  the  data  are  valid,  but 
they  are  much  more  useful  in  qualitative  analysis  of  data  with 
higher  cyclic  tendencies  than  this  material  appears  to  have. 

As  the  value  of  the  mean  of  the  residual  in  the  full  set 
Correction  Factor  test  is  0.00302  (Appendix  2,  page  25)  the  unbiased 
nature  of  the  regression  operation  is  apparent.  The  sum  of  squares 
of  the  residual  is  used  later  in  covariance  analyses  (Appendix  4) 
to  compare  full  data  regression  with  that  of  classifications.  The 

o 

value  given  here  of  0.90350448  x  10°  (Appendix  2,  page  25)  is 
comparable  with  that  given  in  the  last  step  Of  regression  of 
0.90347647  x  10^  (Appendix  2,  page  24,  line  36).  The  value  calcu¬ 
lated  from  the  residuals  is,  however,  somewhat  more  accurate  than 
that  calculated  in  the  repeated  matrix  inversion  of  stepwise  multi¬ 
ple  regression. 
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The  coefficients  of  skewness  and  kurtosis  are  calculated 
from  the  third  and  fourth  moments  of  the  distribution  of  the 
residuals.  If  the  distribution  were  normal  they  should  not  be 
significantly  different  from  zero  and  3.0.  Testing  both  these 
statistics  with  their  given  standard  deviation  it  is  seen  that  the 
distribution  is  positively  skewed  and  somewhat  less  peaked  (kurtotic) 
than  that  which  can  be  ascribed  to  sampling  error. 

In  general  then  we  find  from  test  of  the  residual  that  un¬ 
qualified  acceptance  of  normality  and  independence  is  not  valid. 

If  the  periods  of  measurement  had  been  random  rather  than  fixed 
and  the  number  of  workmen  more  varied  it  is  quite  possible  randomly 
sampled  observations  would  have  produced  more  exactness  in  the 
results.  One  must,  however,  make  the  best  of  what  is  given. 

The  effects  of  non-normality  on  the  results  of  analysis  of 
variance  and  other  related  tests  have  been  investigated  by  a  number 
of  authors  (16,  17). 

It  appears  that,  with  the  degree  of  non-normality  present  in 
these  results,  the  precision  of  confidence  limits  are  decreased 
slightly  but  not  below  usable  levels.  When  a  statistic  has  a  con¬ 
fidence  interval  of  95  per  cent  this  interval  actually  may  be  as 
high  as  the  97  per  cent  or  as  low  as  the  93  per  cent  limits.  Further 
examination  of  particular  parameters  and  their  associated  use  will 
show  that  this  level  of  precision  is  sufficient  for  the  problem 


at  hand. 
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(3)  Multiple  Regression  -  Labour  Content  Estimation 

Model  I 

Model  I  is  characterized  as  follows: 

Y  =  B0  Xc  +  Bx  Xx  +  .  .  .  +  Bk  Xk 

Where  Y  =  Total  contract  hours  per  contract  per  month 
Xo  =  1 

,  X2  .  .  .  X^  =  The  values  of  the  "k"  different 

types  of  production  measured  in 
each  contract 

B-p  B2  .  .  .  =  The  contract  hours  per  unit  of 

production  measured 

The  results  from  the  full  scale  test  of  this  model  using  all 
the  data  (1275  observations)  are  given  in  Appendix  1.  This  compu¬ 
tation  was  a  stepwise  regression  analysis,  as  outlined  in  Theory, 
on  variables  1,2,3,...  29,  30,  32,  and  33.  (Specific  descrip¬ 
tions  of  these  job  types  are  given  in  Table  1) .  Production  unit  31 
was  not  included  as  only  one  value  for  this  unit  was  present  in  the 
observations . 

In  selecting  the  regression  step  (Step  31)  which  produced  the 
"best"  equation  the  following  statistics  were  used  as  selection 
criteria^ 

1.  The  "F"  ratio  of  the  mean  square  for  regression  and  the 

mean  square  for  error  of  the  analysis  at  all  stages  of  regression 
was  highly  significant.  This  ratio  and  the  probability  that  a  value 
as  great  as  this  could  occur  by  chance  are  given  in  columns  4  and  5 
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of  the  analysis  of  variance  following  each  regression  step.  Examin¬ 
ation  of  Program  3  (Appendix  1)  shows  that  "F"  ranges  from  1162.8 
at  step  1  to  210.76  at  step  32.  The  probability  function  of  this 
value  however  is  always  0.00000.  One  can  therefore  assert,  with 
almost  complete  certainty,  that  the  regression  function  at  all  steps 
is  valid  by  analysis  of  variance  test. 

2.  The  quantity  "tn,"  being  the  ratio  of  the  difference 
of  sums  of  squares  for  regression  at  step  ,?n"  and  step  (n-1)  and 
the  mean  square  for  error  at  stage  "n,"  measures  the  contribution 

of  the  variable  Xn  to  the  regression  after  variables  X-^,  X2 ,  ...Xn-1, 

u>*  S 

have  been  added,  wet**  used  to  decide  which  step  of  regression  to  be 
used  as  a  working  equation.  This  value  (tn)  is  distributed  as  "t" 
but  may  be  considered  as  an  exact  normal  distribution.  Any  values 
therefore  of  11  tn”  greater  than  1.96  indicate  that  the  addition  of 
variable  Xn  has  a  significance  greater  than  95  per  cent.  In  Program 
3  we  see  that  "tn"  (line  4^ column  5jOf  the  analysis  of  variance 
tables)  decreases  suddenly  to  a  non-significant  value  of  0.08.  Step 
32  of  the  regression  analysis  is  therefore  disregarded  and  Step  31 
accepted  as  the  working  equation. 

3.  A  second  statistic  which  is  related  to  "tn,"  the  sum 
of  the  percentage  of  variance  reduction^ also  indicates  that  Step  31 
should  be  the  working  equation.  The  values  for  Steps  29  through  32 
are  84.38,  84.42,  84.45  and  84.45.  As  Step  32  of  regression  does  not 
reduce  the  variance  of  the  dependent  variable  below  that  of  Step  31 
the  use  of  this  added  dimension  is  superfluous. 
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4.  Examination  of  the  coefficients  of  the  chosen  regres¬ 
sion  step  (No.  31,  Appendix  1,  page  23)  in  detail  now  was  done. 

The  statistics  used  were  the  "ta(3i)U  values*  This  value,  the 
ratio  of  the  estimated  coefficient  divided  by  its  standard  deviation, 
is  distributed  as  ut"  with  N-n-1  degrees  of  freedom  and  may  be  used 
to  test  a  coefficient  for  a  significant  difference  from  zero.  Again, 
with  over  1200  degrees  of  freedom^ values  less  than  1.96  indicate  that 
there  is  more  than  a  1  in  20  chance  that  the  true  coefficient  could 
be  zero.  Of  the  31  ta^^  values  given  (Appendix  1,  page  19'j,  it 
was  seen  that  only  those  for  variable  32  (Miscellaneous  Timber  -  $) 
and  29  (Miscellaneous  Development)  are  below  this  value.  Consider¬ 
ing  that  these  ta(3i)  values  (1.75  and  1.50)  were  only  slightly  below 
the  95  percent  level  the  retention  of  these  coefficients  seemed 
reasonable . 


5.  A  final  examination  of  the  31  coefficients  showed  that 
all  estimated  values  were  positive.  There  are  no  restrictions  in 
the  model  or  the  computational  procedure  against  negative  coefficients 
but  the  practical  conception  of  a  negative  time  per  production  unit 
would  be  most  difficult  in  any  mine  contract  system.  The  results 
are  thus  realistic  from  this  point  of  view. 

From  the  above  criteria  then  the  following  equation  is  selected 
as  ’’best"  under  the  Labour  Content  Model: 
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Y  =  92.35  +  1.90  Xi  +  1.78  X2  +  1.60  X3 
+  1.65  X4  +  2.43  X5  +  1.50  X6  +  0.149  X? 

+  0.114  X8  +  0.0803  X9  +  0.0608  X10 
+  0.135  XU  +  0.0866  X12  +  0*0822  X13 
+  0.0559  X14  +  0.0271  X15  +  0.0165  X16 
+  2.90  Xly  +  3.97  Xlg  +  0.0911  X19 
+  0.0376  X20  +  0.923  X21  +  0.235  X22 
+  0.221  X23  +  0.741  X24  +  0.550  X25 
+  0.130  X26  +  0.0251  X27  +  0.0679  X28 
+  0.828  X29  +  0.221  X3Q  +  0.161  X32  (7) 

The  dimension  of  Y  is  contract  hours  and  the  coefficients 
are  in  dimensions  of  contract  hours  per  production  unit.  The  values 
are  given  to  three  significant  figures  as  that  is  the  normal  measure-- 
ment  significance  of  production  units.  The  notation  Y  is  given  rather 
than  Y  to  indicate  that  this  is  a  predicted  value  for  total  contract 
hours  per  contract  rather  than  an  exact  mathematical  value. 

To  compare  the  coefficients  of  the  above  equation  with  the 
linear  equation  represented  by  the  contract  system  rates,  which  are 
in  dollars  per  unit  of  production,  one  must  convert  equation  (7)  to 
equivalent  dimensions.  A  second  factor  that  must  be  resolved,  be¬ 
fore  comparisons  can  be  made,  is  that  the  constant  of  92.35  hours 
in  the  regression  equation  is  not  represented  in  the  equation  of  the 
contract  system.  (Constant  non-productive  time  is  certainly  considered 
when  constructing  contract  rates  from  time  study  methods  but  in  the 
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given  case,  rates  have  been  adjusted  upward  to  include  this  factor). 

With  the  above  in  mind  the  following  calculations  were  made 
on  the  values  of  equation  7  and  its  associated  ta^3^  values. 

(a)  Weighting  factor  to  adjust  net  labour  content 
estimates  to  gross  or  full  time  values  — 


Let 

x34  = 

mean  contract  hours  per  contract 

Bo  * 

estimated  constant  contract  hours  per  contract 

Ki  ■ 

constant  for  multiplicative  adjustment  of 
labour  content  estimates  (non-dimensional) 

Then 

ii 

rH 

!*: 

x34 

^34  "  Bo 

And  for 

Program  3 

Ki  ■ 

404 .2  =  1  30 

(404.2  92.35) 

(b)  Conversions  of  the  labour  content  estimates  of  hours 
per  unit  of  production  to  comparable  contract  rates  of 
dollars  per  unit  of  production  - 

Let  S(x|.)  =  the  sum  of  all  dollar  credits  paid  to 

1~*  all  contracts 

Let  S(Y^)  =  sum  of  all  contract  hours  worked 

(i  =  1,  1275,  j  =  1,  2-. ..29,  30,  32) 

s(x2 .) 

Then  Ko  =  weighted  mean  credit  rate  =  _ _U — 

S(Yi) 

Now^ as  means  are  summations  divided  by  N,  the  ratio  of  two 
means  is  equal  to  the  ratio  of  the  summations  if  the  means 
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are  from  the  same  set.  We  therefore  calculate  K2  from 
the  sum  of  the  means  of  the  variates  of  Program  60  to 
which  dollar  conversion  factors  have  been  applied  and 
the  mean  value  of  in  Program  3  giving: 

k2  =  =  $3.58  per  hour 

(This  is  the  average  earnings  per  hour  of  all  men  on 
contract  for  the  period  February  1962  to  November  1964), 

The  two  factors  K  and  K£  can  now  be  combined  to  give  an 
adjusted  earnings  rate  which  when  multiplied  by  the 
labour  content  estimated  for  a  production  unit  gives 
a  rate  in  dollars  per  hour  comparable  to  the  contract  rate. 

Then  Rc  (the  calculated  contract  rate) 

=  (K2)  .  (Kj)  .  (Bj)  .  =  (4.65) (Bj)  dollars/unit 


(c)  The  confidence  limits  of  the  estimates  must  also  be 
adjusted.  The  fractional  method  using  the  ta  statistic 
is  the  easiest  approach. 


As  the  95  per  cent  confidence  limits  for  B.  are  defined 

•Jr 

by  Bj  ±  (1.96)  (Sa), where  (a  =  j)  and  Sa  -  standard  error 

of  estimate  or  standard  deviation  of  B^ ,  we  may  express  this 
confidence  limit  as  a  fraction  of  B^  by  K3  = 


(1.96)  (Sa)  _  (1.96) 
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As  an  example: 

K3(i3)  (for  BU)  -  -  0.109 


Now  (Program  3  Step  31)  =  0,0822 

Then  Rc(13)  =  (4.85)  (0.0822)  =  0.399  $/ton. 


The  95  per  cent  confidence  limits  for  R 


c  (13) 


are 


0.399  +  (0.109)  (0.399)  =  0.399  +  0.043  $/ton. 


By  95  per  cent  confidence  limits  it  is  implied  that  if  we 
compare  Rcq3^  with  the  contract  rate  for  Unit  13  and  if  the 
contract  rate  for  Unit  13  is  greater  than  0.442  dollars  per 
ton  or  less  than  0.356  dollars  per  ton  then  there  is  a 
20  to  1  chance  that  the  contract  rate  is  biased. 

The  above  calculations  were  done  for  all  estimates  of  the 
full  scale  Labour  Content  Model  test  and  are  given  in  Table  4. 

Examination  of  Table  4  indicates  that? in  general,  contract  rates 
have  been  equitable  over  the  three  year  period  of  observations.  The 
rates  which  are  significantly  different  from  those  estimated  as 
"best,"  indicated  by  an  *,  comprise  a  relatively  small  part  of  the 
total  time.  One  also  must  appreciate  that  the  analysis  has  been  made 
on  data  from  a  relatively  long  time  period  which  may  include  some  out 
of  date  conditions.  Before  any  change  of  rates  is  made  on  the  basis 
of  this  analysis  one  should  determine  if  the  apparent  inequalities  of 
rates  are  operative  at  the  present  in  all  methods  of  mining. 
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COMPARISON  OF  LABOUR  CONTENT  ESTIMATES  AND  CONTRACT  RATES 


Variable 

No.  Description 

Rc($/unit) 
Rate  Calc. 

from 

M.  R.  A. 

95%  Conf. 
limits  of 
Rc  + 
$/unit 

Contract 

Rates 

Given 

$/unit 

Unit  of 
Measure 

1. 

Credit  time 

8.85 

1.18 

1.56** 

hours 

2. 

Subdrifting  0*  -  50* 

8.36 

1.10 

6.90* 

feet 

3. 

Subdrifting  50*  -  100' 

7.45 

1.18 

7.40 

feet 

4. 

Subdrifting  100*  -  150' 

7.67 

1.76 

7.30 

feet 

5. 

Subdrifting  +  1508 

11.30 

4.41 

8.40 

feet 

6 . 

Stope  Raising 

6.98 

2.09 

6.40 

feet 

7. 

Breaking  4‘  -  6’  widths 

0.70 

0.24 

0.67 

tons 

8. 

Breaking  6*  -  10 1  widths 

0.53 

0.06 

0.48 

tons 

9. 

Breaking  10*  -  20'  widths 

0.37 

0.05 

0.42 

tons 

10. 

Breaking  +  20*  widths 

0.28 

0.05 

0.39** 

tons 

11. 

Slushing  5  HP  air 

0.63 

0.19 

0.60 

tons 

12. 

Slushing  10  HP  air 

0.40 

0.60 

0.50** 

tons 

13. 

Electric  Slushing 

0.38 

0.04 

0.40 

tons 

14. 

Second  Slush 

0.26 

0.08 

0.20 

tons 

15. 

Slusher  Floor  0*  -  7' 

0.13 

0.04 

0.12 

sq.ft. 

16. 

Slusher  Floor  +  7s 

0.08 

0.02 

0.10 

sq.ft. 

17. 

Cribbed  Millhole 

13.50 

3.14 

15.00 

feet 

18. 

Stulled  Millhole 

18.45 

4.03 

13.50* 

feet 

19. 

Gob  Fence 

0.42 

0.10 

0.30* 

sq.ft. 

20. 

Hangingwall  Fence 

0.17 

0.07 

0 . 30** 

sq.ft. 

21. 

Test  holes 

4.24 

3.46 

0.75* 

each 

22. 

Rock  Bolts  W/Plate 

1.09 

0.70 

2.50** 

each 

23. 

Rock  Bolts  No  Plate 

1.03 

0.98 

1.25 

each 

24. 

Posts 

3.44 

2.14 

2.50 

each 

25. 

Cribs 

2.56 

2.54 

3.00 

feet 

26. 

Bulkhe ad 

0.60 

0.32 

0.40 

sq.ft. 

27. 

Tramming 

0.12 

0.05 

0.15 

tons 

28. 

Muck  &  Tram 

0.31 

0.08 

0.38 

tons 

29. 

Mi  sc.  Development 

3.85 

5.05 

1.00 

dollars 

30. 

Mi  sc.  Breaking 

1.03 

0.59 

1.00 

dollars 

31. 

Misc.  Timber 

0.75 

0.84 

1.00 

dollars 

NOTE:  *  indicates  rates  differing  beyond  the  95%  limits 
**  indicates  those  beyond  the  99%  limits 
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An  answer  to  the  above  is  attempted  in  subsequent  regression 
analyses  on  subsets  of  the  data. 

(4)  Multiple  Regression  -  Correction  Factor  Estimation 

Model  II 

The  results  of  the  full  scale  test  (Program  60)  of  this  model 
on  all  the  data  (1275  observations)  are  given  in  Appendix  2.  The 
independent  variables  in  this  model  are  derived  from  the  same  raw 
data  as  used  in  the  former  labour  content  estimation  model  by  multi¬ 
plication  with  the  contract  rates.  The  factor  used  to  convert  the 
observed  contract  hours  to  dollars  was  $3 . 6025/hour .  This  factor 
was  estimated  from  a  number  of  trial  runs  and  represents  a  weighted 
mean  earnings  rate  per  observation  (contract  calculation) .  This 
rate  is  somewhat  greater  than  the  mean  of  bonus  rates  (see  variable 
No.  41  Appendix  1,  page  1)  plus  the  standard  contract  wage  debit 
rate  of  $1.56  per  hour.  Some  of  this  descrepancy  can  be  attributed 
to  penalty  charges  but  it  is  quite  possible  that  contract  calcula¬ 
tions  with  high  total  hours  have  tended  toward  higher  bonus  rates. 
This  should  be  expected  as  the  larger  more  stable  contracts  would 
tend  to  attract  the  senior  and  more  skilled  workmen. 

Examination  of  the  output  of  Program  60  shows  that  while  the 
coefficients  estimated  differ  from  those  of  Program  3  the  statistics 
used  to  evaluate  the  regression  and  the  parameters  within  the  regres¬ 
sion  are  exactly  the  same  as  those  of  Program  3.  The  reason  for  this 
apparent  incongruity  is  that  statistics  such  as  means  and  regression 
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coefficients  are  absolute  values  while  test  statistics  such  as  the 
"F"  ratio,  percent  variance  reduction,  "tn"  and  "t  "  are  calculated 
from  ratios  in  which  multiplicative  factors,  applied  to  the  data, 
cancel . 

In  order  to  demonstrate  that  the  two  methods  of  approach  are 
equivalent,  a  method  was  developed  to  produce  contract  rate  cor¬ 
rection  factors  from  the  estimates  of  Program  3.  The  assumption 
behind  this  method  is  that, if  given  a  set  of  just  contract  rates, 
the  percent  of  time  spent  in  producing  any  given  measured  units 
should  be  equal  to  the  percent  of  money  paid  for  this  production. 
The  calculation  method  follows  for  the  results  which  are  tabulated 
in  Table  5. 

(i)  Distribution  of  dollar  credits  for  each  unit  "j"  of 
production  as  a  percent  of  total  contract  credits: 


(100) 


s  X 


where  is  the  mean  of  the  "i"  transformed  variable 

in  Program  60, and  S  indicates  summation  over  all  the 
means  in  the  same  set. 

(ii)  Distribution  of  time  for  each  unit  "j"  of  production 
as  a  percent  of  total  contract  time: 
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— r —  (ki)  (100) 

S(xp 

Where  X^  is  the  mean  value  for  unit  "  j",  and  HS"  indicates 
summation  of  over  all  means, and  is  as  defined  for  Table  4 
calculations . 

(iii)  Correction  factor  calculation  from  Program  3 
=  7»  mean  time  +  %  mean  credits 

(iv)  Correction  factors  calculated  from  Program  60 


Where  B.  is  the  regression  coefficient  from  Program  60, 

J  _i 

Step  31,  and  is  the  mean  value  of  the  total  dollar 

credits,  the  dependent  variable  of  Program  60. 

(v)  The  percentage  confidence  limits  shown  in  the  last 
column  of  Table  5  are  derived  directly  from  the  fractional 
confidence  limits . 

(5)  Comparison  of  Results  from  Model  I  and  Model  II 
Examination  of  correction  factors  calculated  from  the  two 
regression  models  (Table  5,  columns  3  and  4)  shows  that  they  pro¬ 
duce  almost  identical  results.  Except  for  rounding  errors  the 
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TABLE  5 

COMPARISON  OF  LABOUR  CONTENT  ESTIMATION 
AND  CORRECTION  FACTOR  ESTIMATION 


Variable 

No.  Description 

%  Total 
Credits 

% 

Total 

Time 

Calc. 
Corr . 
Factor 
Prog.  3 

Corr. 
Factor 
Prog.  60 

95% 
Limits 
as  7. 

1. 

Credit  Time 

0.6 

3.7 

6.17 

5.70 

13.3 

2. 

Subdrifting  0-50' 

5.8 

6.9 

1.19 

1.21 

13.1 

3. 

Subdrifting  50'  -  100* 

4. 6 

4.6 

1.00 

1.01 

15.9 

4. 

Subdrifting  100'  -  150' 

1.7 

1.7 

1.00 

0.97 

23.0 

5. 

Subdrifting  +  150' 

0.4 

0.6 

1.50 

1.35 

39.4 

6. 

Stope  Raise 

1.0 

1.1 

1.10 

1.09 

30.0 

7. 

Breaking  4*  -  6' 

1.2 

1.2 

1.00 

1.04 

33.7 

8. 

Breaking  6'  -  10' 

16.9 

18.3 

1.08 

1.11 

11.4 

9. 

Breaking  10'  -  20' 

7.1 

6.3 

0.83 

0.89 

14.7 

10. 

Breaking  +  20' 

6.2 

4.5 

0.73 

0.73 

17.3 

11. 

Slushing  5  HP  air 

0.8 

0.8 

1.00 

1.05 

31.2 

12. 

Slushing  10  HP  air 

3.0 

2.4 

0.80 

0.81 

22.6 

13. 

Electric  Slushing 

24.9 

23.6 

0.95 

0.96 

10.9 

14. 

Second  Slushing 

1.4 

1.8 

1.29 

1.30 

30.5 

15. 

Slusher  Floor  O'  -  7' 

1.8 

1.9 

1.06 

1.05 

27.3 

16. 

Slusher  Floor  4-  7' 

2.5 

2.0 

0.80 

0.77 

29.6 

17. 

Cribbed  Millhole 

3.9 

3.5 

0.90 

0.90 

23.2 

18. 

Stulled  Millhole 

1.7 

2.3 

1.36 

1.37 

21.8 

19. 

Gob  Fence 

1.7 

2.4 

1.41 

1.42 

25.3 

20. 

Hangingwall  Fence 

2.6 

1.5 

0.58 

0.59 

42.5 

21. 

Test  Holes 

0.1 

0.6 

6.00 

5.75 

82.0 

22. 

Rock  Bolts/Plate 

1.4 

0.8 

0.57 

0.56 

64.1 

23. 

Rock  Bolts  No  Plate 

2.2 

1.4 

0.60 

0.64 

95.5 

24. 

Posts 

0.5 

0.8 

1.60 

1.38 

62.2 

25. 

Cribs 

0.3 

0.3 

1.00 

0.86 

39.6 

26. 

Bulkhead 

0.4 

0.7 

1.75 

1.51 

54.2 

27. 

Tramming 

2.6 

2.0 

0.77 

0.78 

39.1 

28. 

Muck  &  Tram 

2.1 

1.8 

0.86 

0.83 

25.5 

29. 

Misc.  Development 

0.0 

0.0 

-- 

3.87 

132.0 

30. 

Misc.  Breaking 

0.3 

0.3 

1.00 

1.03 

52.3 

31. 

Misc.  Timber 

0.3 

0.2 

0.75 

0.75 

112.0 
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factors  are  the  same.and  as  the  t  ,  N  values  are  identical  the 
same  nine  production  rates  "1,  2,  10,  12,  18,  19,  20,  21  and  22” 
appear  inappropriate.  It  can  therefore  be  concluded  that 
regression  analysis  Model  II  is  a  valid  method  to  use  for  analysis 
of  this  data.  The  advantages  of  this  method  over  Model  I  are  apparent 
when  the  calculations  required  to  produce  Tables  4  and  5  are  considered, 
Whereas  Column  3  of  Table  5  requires  several  steps  of  multiplication 
and  summation.  Column  4  of  Table  5  consists  of  only  one  multiplica^ 
tive  step  beyond  that  given  in  the  computer  output.  As  the  raw 
estimates  in  Program  3  are  in  several  different  units,  direct  inter¬ 
comparison  of  estimates  is  impossible.  In  Program  60  all  raw  estimates 
are  in  the  same  units  and  one  can  easily  compare  pairs  or  larger  sets 
of  estimates  directly.  A  number  of  multiple  range  tests  are  avail¬ 
able  (18)  for  this  type  of  comparison. 


(C)  Further  Test  to  Establish  the  Effectiveness  and 

Reliability  of  Model  II 

(1)  General  Requirements  for  an  Effective  Contract 

Administration  Analysis  Tool 

It  has  been  demonstrated  in  the  foregoing  that  an  analysis 
method  is  available  which  appears  to  give  a  reasonable  and  unbiased 
evaluation  of  the  past  three  years  contract  data.  It  must  now  be 
shown  that  this  method  would  be  applicable  to  present  and  future 
administration  and  control  of  this  contract  system. 
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The  major  question  that  arises  in  the  consideration  of  any 
contract  system  is  the  equity  of  the  rates.  The  level  of  incentive 
pay,  the  classification  of  units,  the  type  and  timing  of  measure¬ 
ment  and  the  methods  of  calculation  are  generally  fixed  over  relatively 
long  periods  of  time.  The  contract  systems  are  usually  constructed 
so  that  changes  in  work  methods  and  conditions  can  be  reflected  in 
rate  changes .  Obtaining  an  unbiased  and  controllable  measure  of 
what  these  changes  should  be  has  always  been  a  problem.  The  coef¬ 
ficients  from  the  Correction  Factor  multiple  regression  give  unbiased 
estimates  of  a  multiplicative  factor  for  rate  correction.  More  im¬ 
portant,  these  estimates  have  measures  of  reliability,  which  is 
more  than  can  be  said  for  many  time  study  methods. 

A  further  question  of  application  of  these  estimates  is  de¬ 
pendent  upon  stability.  The  data  used  in  Programs  3  and  60  covered 
slightly  less  than  three  years  operation.  The  estimates  can  be 
thought  of  as  mean  corrective  factors  that  could  have  been  applied 
over  this  period  and  should  be  applied  for  further  work  if  conditions 
remain  the  same. 

No  reasonable  administrator  attempts  to  foretell  the  future 
exactly.  The  best  he  can  hope  for  is  to  note  from  the  immediate 
past  that  conditions  are  changing  and  then  try  to  discover  what  is 
causing  the  change  by  use  of  judgement  and  intuition  distilled 
from  experience.  With  the  above  in  mind  one  can  appreciate  the 
value  of  time  series  classification  for  analysis  using  Model  II. 
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If  conditions  and  work  methods,  as  analysed  by  Model  11^ 
appear  sufficiently  stable  in  time,  for  control  purposes  one 
still  must  investigate  other  areas  of  possible  instability.  It  is 
possible  that  the  contract  system’s  assumption  that  work  conditions 
are  relatively  equivalent  could  be  in  error.  A  number  of  regres¬ 
sions  with  the  data  classified  according  to  differing  work  methods 
should  allow  a  test  of  this  assumption. 

A  third  factor,  relative  incentive  rate,  may  also  be  opera¬ 
tive  in  work  rate  stability.  Are  the  work  rates  of  workmen  who 
consistently  earn  better  than  average  incentive  rates  relatively 
different  than  those  of  low  and  average  incentive?  Rate  classes 
should  answer  this  question. 

The  same  model  as  used  in  Program  60  is  used  to  test  the 
above  three  questions  of  stability.  The  methods  of  classifica¬ 
tion  have  been  given  in  Table  2.  (Appendix  3, pages  6  and  7). 

The  computer  outputs  for  these  calculations  were  bound  and  are  on 
file  at  the  Department  of  Mining  and  Metallurgy,  University  of 
Alberta,  Edmonton,  Alberta  (19,  20,  21).  Covariance  analyses  to 
determine  the  validity  of  the  classifications  are  given  in 
Appendix  4. 

Examination  of  the  covariance  analyses  indicates  quite  defi¬ 
nitely  that  all  three  classifications  are  valid.  The  large  "F” 
ratios  indicate  that  the  reduction  of  variance  by  classification  is 
very  much  greater  than  that  which  could  be  attributed  to  chance. 
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(2)  Organization  of  Results  of  Class  Regressions 

for  Comparative  Analysis 

The  thorough  analysis  of  large  quantities  of  data  is 
often  a  problem  to  the  experimenter  and  the  administrator.  Tabular 
lists  of  figures  as  produced  directly  by  computer  are  only  compre¬ 
hensible  to  a  degree.  With  a  computer's  facility  for  producing 
these  tables  one  is  easily  led  to  produce  such  a  confusion  of  figures 
and  details  that  the  solution  is  more  confusing  than  the  problem. 
Examination  of  the  computer  results  from  regressions  on  a  few  classi¬ 
fications  indicated  that  the  likelihood  of  this  type  of  confusion 
was  high.  As  a  method  of  condensing  results  from  a  large  number  of 
regressions  was  required,  it  appeared  that  the  Autoplot  system  (13) 
offered  the  best  solution. 

A  FORTRAN  program  was  written  to  accept  the  parameters  of 
selected  sets  of  regressions  and  present  them  to  the  Autoplot 
system  in  a  form  similar  to  industrial  control  chart  methods.  The 
coefficient  estimates,  (central  line  on  the  plots)  were  pro-rated 
to  a  full  time  basis  in  the  same  manner  as  the  coefficients  in 
Table  4.  The  weighting  factors  varied  from  1.1  to  1.3.  However,  the 
confidence  limits  (outer  lines  on  the  plots)  were  calculated  from 
1.96  times  the  unweighted  standard  error  of  estimate  of  the  coef¬ 
ficient  and  hence  do  not  represent  a  95  per  cent  confidence  band 
but  a  band  of  85  per  cent  to  95  per  cent  confidence. 
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The  step  of  regression  chosen  as  a  working  equation  for  plot¬ 
ting  was  selected  using  the  t  statistic  as  the  major  selection 
criterion.  As  the  classification  regression  sets  are  made  up  of 
smaller  numbers  of  observations  than  the  full  set,  the  selected 
number  of  significant  coefficients  is  always  less  than  that  of  the 
full  set.  It  was  found  that  t^  values  generally  decreased  sharply 
below  1.96  at  the  24th  to  28th  step  of  regression. 

The  combination  of  all  valid  regression  coefficients  in  a 
number  of  plots  (one  for  each  production  unit)  indicated  that  in  a 
number  of  cases  estimates  for  the  coefficient  did  not  appear  in  all 
regressions.  It  was  decided  therefore  to  retain  only  those  plots 
in  which  at  least  75  per  cent  of  all  classes  were  represented.  The 
regression  data  from  which  the  plots  are  made  is  on  file  as  refer¬ 
ence  in  the  Department  of  Mining  and  Metallurgy,  (19,  20,  21). 

The  scales  on  these  Autoplot  graphs  may  at  first  glance  be 
confusing  as  they  are  given  in  exponential  form.  Note  that  on  the 
left  of  the  axes  one  or  two  numbers  are  given  preceded  by  a  sign 
and  a  letter  "E."  This  represents  the  power  of  10  by  which  a 
particular  scale  is  multiplied  to  produce  its  true  value.  For 
example  Figure  50  has  E-01  above  the  horizontal  and  left  of  the 
vertical  axis  and  also  below  the  horizontal  and  left  of  the  vertical 
axis.  These  values  indicate  that  both  scales  should  be  multiplied 
by  10"1. 
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In  consideration  of  these  plots  one  must  appreciate  that 
the  only  true  data  points  are  at  the  classification  points  and 
all  other  points  have  been  produced  by  straight  line  extrapolation 
between  contiguous  classifications.  In  the  plots  of  time  classifi¬ 
cation  and  in  the  incentive  rate  classification  plots  this  extra¬ 
polation  allows  one  to  visualize  trends.  In  the  plots  of  mining 
method  classification  however  one  should  regard  the  extrapolation 
as  extraneous  and  visualize  the  plots  much  the  same  as  a  bar  chart. 

In  interpretation  of  the  plots,  in  the  following  sections, 
reference  will  be  made  to  productivity  increases  or  decreases,  as 
shown  by  the  fluctuations  of  the  correction  factor.  When  the  cor¬ 
rection  factor  for  a  contract  rate  is  above  1.00  productivity  can 
be  said  to  have  decreased  below  that  implied  by  the  given  rate. 

Similarly  when  the  factor  falls  below  1.00  one  can  say  productivity 
has  increased. 

(3)  Analysis  of  Individual  Parameters  within  Classifications 

by  Means  of  Control  Chart  Methods 

(a)  Classification  on  the  Basis  of  Time  Intervals 

A  six  month  moving  average  calculated  at  two  month  in¬ 
tervals  is  used.  This  method  was  chosen  to  obtain  a  balance  between 
the  poor  estimation  that  could  result  from  the  fewer  observations  in 
a  shorter  period  and  the  lack  of  variation  in  longer  periods.  If 
this  type  of  analysis  were  to  be  applied  as  a  control  method  several 
test  analyses  should  be  done  with  varying  length  periods  to  ascertain 
the  method  most  applicable  to  the  available  data  and  methods  of  control. 
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The  plots  produced  by  this  method  of  analysis  generally 
show  good  evidence  of  regularity.  The  overlapping  method  of 
classification  would  produce  some  smoothing  but  it  appears  that 
an  overriding  evidence  of  stability  is  present.  Another  factor 
which  reinforces  the  conclusion  of  stability  is  the  small  varia¬ 
tion  in  the  confidence  band  widths.  With  unstable  estimates  one 
would  also  expect  the  standard  deviations  of  the  estimates  to 
fluctuate  wildly. 

This  type  of  analysis,  time  series  classification,  is 
meant  to  simulate  the  method  a  contract  administrator  would  employ 
on  a  routine  monthly  basis.  In  a  real  situation  the  administrator 
would  be  in  intimate  contact  with  the  system  and  should  have  much 
pertinent  information  for  the  explanation  of  fluctuations.  The 
author  has  had  some  association  with  the  system  under  analysis 
and  will  put  forward  some  causes  of  fluctuations  to  illustrate 
the  method's  use  but  one  must  appreciate  that  fuller  use  of  the 
method  could  be  made  by  those  who  were  directly  concerned  with 
the  real  system. 

(i)  Mean  Bonus  Rate  and  Mean  Contract  Hours 
vs  Time  Figures  16  and  17 

These  mean  values  were  plotted  with  the 

thought  that  they  might  help  explain  some  of  the 

fluctuations  in  the  regression  estimates.  If  a 

decrease  in  the  work  rate,  for  a  production  unit, 

accompanies  an  increase  in  bonus  rate,  a  reduction 


in  contract  rate  would  be  somewhat  validated. 
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The  general  appearance  of  each  plot  suggests  that 
productivity  has  shown  a  general  improvement  over  the 
three  year  period.  The  number  of  contract  hours  per 
contract  also  shows  evidence  of  seasonal  trend.  This  could 
be  a  reflection  of  variation  in  work  force  due  to  seasonal 
holidays . 

(ii)  Percent  Variance  Reduction,  Weighting 
Factor  and  Mean  of  Dollar  Credits  vs 
Time  -  Figure  18 

These  three  factors  appear  in  Figure  18  from  top 
to  bottom  as  listed  above.  The  Weighting  Factor  has 
been  multiplied  by  100  and  the  Dollar  Credits  divided 
by  10.  Percent  Variance  Reduction  is  in  true  scale  and 
is  a  measure  of  the  validity  of  the  regressions  through¬ 
out.  The  Weighting  Factor  is  a  measure  of  the  constant 
in  the  regression  and  may  be  thought  of  as  a  measure  of 
the  non-productive  contract  time.  The  Mean  Dollar  Credits 
are  a  direct  measure  of  mean  total  contract  hours  per 
contract, as  a  constant  multiplier  ($3.6025)  was  used 
to  convert  the  value  of  observed  total  hours  to  dollars. 

The  Variance  Reduction  plot  shows  good  stab¬ 
ility  and  indicates  that  all  regression  equations  were 
highly  significant.  The  values  for  this  parameter  be¬ 
tween  90  per  cent  and  79  per  cent  are  quite  comparable 
to  that  (84.5%)  obtained  from  the  complete  set  of  data 
(Program  60,  Appendix  2,  page  23). 
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The  Weighting  Factors  which  were  applied  to 
all  coefficients  in  each  regression  appear  somewhat 
less  stable,  but  except  for  one  point  all  values  are 
between  1.18  and  1.45. 

The  mean  value  of  dollar  credits  give  a  smooth 
curve  with  a  maximum  in  the  last  half  of  1963.  This 
could  indicate  relatively  larger  crews  on  contracts 
in  this  period  or  it  could  be  caused  from  increased 
productivity. 

(iii)  Credit  Time  -  Correction  Factor  vs  Time 

Figure  19 

The  plot  of  the  correction  factor  for  Credit 
Time  is  relatively  smooth.  It  shows  however,  a  wide 
range  of  values  (2.0  -  9.0).  A  probable  cause  for 
this  fluctuation  might  be  attributed  to  the  adminis¬ 
tration  of  this  factor.  Credit  time  was  originally 
meant  to  be  paid  for  types  of  work  for  which  no  con¬ 
tract  rates  had  been  set.  Supervisors  were  to  esti¬ 
mate  the  time  the  workmen  spent  doing  this  work  and 
give  sufficient  credit  hours  at  $1.56  per  hour  to 
return  average  incentive  rates.  Often,  however, 
credit  time  may  have  been  confused  with  non-contract 
time  which  was  to  be  applied  only  when  the  workmen,  due 
causes  beyond  their  control,  were  not  able  to  produce. 
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The  overall  effect  of  this  improper  allow¬ 
ance  may  have  been  to  penalize  the  workmen  for  doing 
non-standard  work.  The  effect,  however,  would  not  be 
great,  as  mean  credit  time  comprises  less  than  1.0  per 
cent  of  mean  Credits.  It  is  also  gratifying  to  note  that 
5.70,  the  factor  estimated  with  all  the  data  (Table  5$ 
Column  4?line  1)  is  very  near  the  mean  of  the  plotted 
values.  From  the  foregoing  one  could  suggest  that, 
rather  than  change  the  rate  for  credit  time,  a  clari¬ 
fication  of  administration  policies  would  materially 
assist  in  correction  of  this  anomaly. 

(iv)  Development  Correction  Factor  vs  Time 
Figures  20  -  24 

Overall  examination  indicates  that  develop¬ 
ment  rates  are  equitable  since  most  confidence  bands 
enclose  1.0  throughout.  Subdrifting  0!  -  50'  from 
the  mill  hole  shows  some  trend  to  higher  productivity 
while  similar  work  50'  -  100'  from  the  disposal  point 
trends  in  the  opposite  direction.  Investigation  of 
these  production  units  could  be  in  order.  Equip¬ 
ment  allocation  and/or  work  methods  are  likely  causes 


of  the  trends . 


(v)  Ore  Breaking  Correction  Factors  vs  Time 

Figures  25  -  28 

The  non-continuity  of  the  factors  for  4'  -  6' 
widths  (Figure  58)  is  to  be  expected  as  very  few  stope 
are  mined  at  this  width  alone.  The  correction  factor 
for  the  major  production  width  of  6'  -  10*  (Figure  26) 
is  continuous  and  has  a  relatively  narrow  confidence 
band.  A  definite  low  in  productivity  late  in  1962 
is  followed  by  a  regular  increase  to  the  end  of  1964. 
The  other  two  categories  of  10'  -  20'  and  +  20'  widths 
also  show  this  increasing  productivity  trend.  An  in¬ 
vestigation  of  work  methods  and  equipment  changes  for 
all  ore  breaking  work  may  be  in  order.  Changes  of 
rates  to  reflect  changes  in  productivity  of  these 
units  could  have  a  marked  effect  on  contract  costs  as 
they  account  for  over  31  per  cent  of  costs. 

If  one  compares  the  time  of  increasing  produc¬ 
tivity  in  these  major  breaking  variables  with  the  time 
of  increase  in  bonus  and  decrease  in  contract  hours 
(Figures  16  -  17)  it  seems  that  work  methods  in  ore 
breaking  have  been  improved  toward  the  point  where 
mine  management  should  be  recovering  some  portion 


of  the  improvement. 
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(vi)  Slushing  Correction  Factors  vs  Time 
Figures  29  -  32 

All  figures  except  No.  31,  electric  slushing, 
show  wide  confidence  bands  that  enclose  the  value  of 
1.0.  Hence  no  change  of  rates  seems  warranted.  Electric 
slushing,  the  major  production  variable  in  the  contract 
system,  (24  per  cent  of  costs)?  shows  (Figure  31)  indica¬ 
tions  of  increasing  productivity. 

Changes  of  rates  for  this  job  could  have  a 
decided  effect  on  contract  costs  and  should  be  investi¬ 
gated  closely. 

(vii)  Timbering  Correction  Factors  -  Figures  33  -•  37 
The  plots  of  the  two  types  of  s lusher  floor 
(Figures  33  -  34)  show  their  proper  position  in  time. 

Work  methods  were  altered  in  1963  and  a  lower  rate  was 
established.  It  appears  that  even  with  the  reduced 
rate  much  of  the  profit  from  improvement  in  the  method 
has  accrued  to  the  workmen. 

Millhole  construction  (Figures  35  -  36)  shows 
that  while  cribbed  millholes  have  been  overpriced, 
stulled  millholes  have  been  underpaid.  Changes  in 
these  rates  would  loom  large  to  the  workmen  but 


would  have  little  effect  on  overall  costs. 


■ 
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The  remaining  factor  plotted  in  this  series, 

Gob  Fence  (Figure  37)?  shows  indications  of  over¬ 
pricing  in  the  past  but  is  in  line  with  rates  at 
the  present  time  (1964) . 

(b)  Classifications  on  the  Basis  of  Methods 

of  Work 

Correction  factors  for  each  unit  of  production  were 
plotted  against  work  method  classification.  The  arrangement  of 
the  work  method  classes  in  the  plots  is  somewhat  toward  increas¬ 
ing  irregularity  but  is  in  no  sense  gradational.  While  the  ir¬ 
regularity  in  ore  bodies  has  dictated  some  pecularities  of  each 
method  many  more  factors  have  been  considered  in  the  selection  of  a 
particular  mining  system.  In  examining  the  plots,  one  should  re¬ 
gard  extrapolation  between  classes  as  a  visual  aid  only  and  not  an 
intimation  of  continuity. 

It  will  be  noted  that  the  confidence  bands  tend  to  widen 
for  D  and  E  classes  of  stopes.  As  the  number  of  samples  in  these 
classes  was  only  about  one  third  of  those  in  the  other  classes,  a 
larger  variation  in  estimates  may  be  expected. 

Some  50  samples  classified  as  F  or  S  (unclassified 
strays)  were  not  plotted  as  this  small  highly  varied  sample  returned 
only  a  few  valid  estimates. 

(i)  Percent  Variance  Reduction,  Weighting 

Factor  and  Mean  Dollar  Credits  vs  Stope 

Classification  -  Figure  38 


These  values  are  arranged  on  the  plot  the  same  as 
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similar  values  were  in  the  time  classification  plot 
(Figure  18) . 

The  percent  variance  reduction  graph  is  regu¬ 
lar  throughout.  The  range  of  75  per  cent  to  90  per 
cent  is  similar  to  that  of  the  time  series  classifi¬ 
cation,  an  indication  of  valid  regression  functions. 

The  plot  of  mean  dollar  credits  indicates 
that  A  and  B  types  have  had  considerably  larger  total 
contract  payments  than  other  classes.  This  could  be 
expected  as  more  regular  working  places  are  more 
amenable  to  double  shift  operation. 

The  Weighting  Factor  is  fairly  constant  around 
1.3  except  for  the  C  method.  The  discrepancy  is 
probably  a  reflection  of  the  lower  percent  variance 
reduction  in  the  C  method.  Overall  we  can  expect 
that  the  confidence  limits  shown  on  the  plots  are  in 
the  80  per  cent  to  90  per  cent  range. 

(ii)  Correction  Factor  for  Credit  Time  vs  Stope 

Classification  -  Figure  39 

This  plot  shows  that  the  assessment  of  credit 
time  also  varies  with  mining  method.  Variations  aside, 
there  is  definite  evidence  that  this  factor  has  been 


undervalued. 


' 
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(iii)  Correction  Factors  for  Development  vs 
Stope  Classifications  -  Figures  40  -  42 

These  three  plots  offer  a  possible  clue  to 

the  recent  decreases  in  sub-drifting  productivity  noted 

in  the  time  series  analysis.  "D"  methods  of  development 

work  definitely  have  taken  more  labour  time  than  normal. 

This  type  of  stoping  is  also  of  recent  origin  and  could 

be  the  cause  of  the  general  decline  in  development 

productivity  (Figures  21  -  22) . 

(iv)  Correction  Factors  for  Ore  Breaking  vs 
Stope  Classifications  -  Figures  43  -  44 

It  appears  the  rates  for  breaking  6*  -  10' 

widths  have  been  equitable  throughout  the  mine.  There 

is  some  evidence  that  B  and  C  methods  of  stoping  have 

easier  breaking  in  larger  widths  than  do  the  other 

methods.  Work  methods  and  materials  here  would  bear 

investigation. 

(v)  Correction  Factors  for  Slushing  vs  Stope 
Classifications  -  Figures  45  -  46 

It  appears  that  A  and  D  types  of  stoping  have 
relatively  low  productivity  in  electric  slushing  while 
C  and  E  types  are  relatively  high.  It  is  possible  that 
examination  of  mine  equipment  would  show  a  higher  per¬ 
centage  of  3  drum  slushers  in  places  using  "C"  methods 
than  in  those  using  A  and  D.  The  E  method  of  stoping 
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enjoys  the  unique  advantage  that  workmen  muck  and 
tram  slushed  ore  directly  and  thus  have  no  delays  due  to 
chute  hang -ups. 

It  is  interesting  to  note  that  while  covariance 
analysis  indicates  that  this  classification  of  data  is 
quite  valid  fewer  continuous  estimates  were  calculated. 
Because  many  of  the  production  units  are  particular  to 
one  or  two  methods  of  mining,  measurements  of  these 
units  are  sparse  or  missing  in  many  of  the  classifica¬ 
tions  and  hence  did  not  enter  the  regression  as  valid 
parameters . 

It  also  would  have  been  interesting  to  re¬ 
classify  the  data  under  some  of  the  other  classifications 
such  as  ore  control,  ground  control  or  breaking  charac¬ 
teristics  but  time  did  not  permit. 

(C)  Classification  on  the  Basis  of  Incentive  Rate 

By  this  classification,  an  attempt  is  made  to  ascertain  whether 
work  rates  vary  more  or  less  directly  with  the  incentive  rates.  Do 
the  workmen,  who  consistently  receive  higher  incentive,  work  at  higher 
productive  rates,  on  all  types  of  work,  or  does  the  rate  of  change  of 
productivity  vary  from  job  type  to  job  type? 

The  plotted  values  were  obtained  from  data  sorted  into  $0. 50/hr. 
bonus  rate  increments.  The  estimate  from  the  full  set  regression  has 
also  been  included  and  is  plotted  at  $1. 84/hr.  The  inclusion  of  these 
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values  has  destroyed  much  of  the  continuity  of  the  plots.  Covari¬ 
ance  analysis  on  this  classification  (Appendix  4,  page  3)  and  per¬ 
cent  variance  reduction  indicate  that  considerably  "better"  regres¬ 
sion  equations  are  produced  by  this  classification.  However,  when 
one  considers  that  incentive  rate  values  are  the  core  of  the  contract 
system  it  is  apparent  that  this  classification  is  an  arbitrary  trun¬ 
cation  of  the  data  which  could  materially  degrade  its  essential 
randomness.  This  degradation  could  have  a  decided  effect  on  para¬ 
meter  estimation  and  nullify  any  qualification  of  "better."  One 
therefore,  considers  the  plots  more  for  continuity  and  trends  rather 
than  absolute  level  of  productivity  factors. 

(i)  Percent  Variance  Reduction  and  Mean  Dollar  Credits 

vs  Incentive  Rate  Classification  -  Figure  47 

The  percent  variance  reduction  (the  upper  line  in  Figure  47) , 
a  measure  of  individual  regression  validity,  is  quite  high  and  is 
probably  caused  by  truncation  of  the  data.  The  mean  total  credits 
(lower  line)  show  the  linear  trend  which  could  be  expected  from 
the  classification  method. 

A  Weighting  Factor  was  used  in  plotting  the  variables  but 
was  not  plotted  itself.  Examination  of  the  results  (21)  showed 
it  was  nearly  constant  at  1.3. 

(ii)  Correction  Factors  for  Development,  Breaking, 

Slushing  and  Timbering  -  Figures  49  -  50 

All  of  the  estimates  in  the  above  categories  showed  no 
strong  evidence  of  departure  from  a  straight  line.  One  can  conclude 
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therefore,  that  the  contract  rates  paid  for  these  production  units 
are  relatively  equitable  throughout  the  incentive  rate  range.  (Two 
samples  of  the  20  plots  in  this  classification  are  given  in  Figures 
49  and  50) . 

How  can  one  reconcile  the  above  evidence  with  former  evi¬ 
dence  of  non-equity  in  the  rate  structure?  It  appears  that  the 
mining  methods  are  such  that  overpriced,  equitable  and  under-priced 
units  of  production  tend  to  be  distributed  evenly  throughout  the 
mine.  The  chances  of  a  contractor  consistently  producing  a  surplus 
of  overpriced  units  hence  making  consistently  higher  bonus  has  been 
quite  unlikely. 

(D)  Multiple  Regression  Analysis  of  Expendable  Materials 

Some  of  the  materials  used  in  stope  production  have  been 
recorded  in  the  contract  records.  It  was  decided  to  test  this 
materials  data  with  the  methods  of  regression  that  were  used  to 
evaluate  the  labour  factors  to  find  if  similar  factors  for  material 
use  could  be  obtained.  As  former  regressions  for  labour  content  on 
subsets  have  been  shown  to  be  significant  within  classes,  regres¬ 
sions  for  material  have  been  done  within  classes  only.  Moreover, 
since  materials  are  herein  considered  subsidiary  to  labour  content, 
only  the  mining  method  classification  was  considered  for  illustra¬ 
tion  of  the  method. 

The  five  materials  factors  recorded  in  the  data  were: 

1.  X36  -  Blast  holes  drilled  (fuse  count) 
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2.  X37  -  Reamed  holes 

3.  X38  -  Pounds  of  explosive  as  reported  by  workmen 

4.  X39  -  Pounds  of  high  explosive  as  calculated  by 

supervisor 

5.  X40  -  Pounds  of  Amex  as  calculated  by  supervisor 
(Holes  drilled  for  blasting  and  reamed  holes  are  here  considered 
as  materials  because  there  is  a  direct  relationship  between  the 
number  of  holes  drilled  and  the  amount  of  material  used  such  as  drill 
steel  and  bits) . 

All  the  above  five  variables  were  used  as  dependent  in  separate 
regression  analyses  against  all  measured  production  units  of  jobs  in 
which  they  were  used.  The  computer  outputs  of  five  regressions  in 
all  classes  of  mining  methods  are  given  in  (22) .  The  results  which 
appeared  valid  and  continuous  were  put  in  the  Autoplot  program  in  a 
manner  similar  to  the  labour  factors. 

It  was  found  that  only  two  of  the  dependent  variables  were 
predicted  regularly  and  validly  by  regression  methods.  These  were: 

1.  Blast  holes  drilled 

2.  Pounds  of  explosive  reported 

If  the  two  types  of  explosive  as  calculated  by  supervision 
could  have  been  additive ly  combined  they  might  have  produced  a  valid 
regression.  As  it  was,  it  appears  that  these  two  values  separately 
are  of  little  use.  A  regression  with  these  two  as  independent  vari¬ 
ables  against  holes  drilled  as  a  dependent  variable  could  give  one  an 
estimation  of  the  loading  factor  in  average  pounds  of  explosive  per 


foot  of  blast  hole  drilled. 
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Regression  using  the  number  of  reamed  holes  as  a  dependent 
variable  was  disappointing.  The  use  of  reaming  in  stope  development 
work  is  apparently  much  more  random  than  standard  practice  policies 
would  indicate.  On  the  other  hand  standard  mining  practices  followed 
by  substandard  reporting  policies  could  have  been  a  cause  of  the  above 
anomaly. 

(1)  Holes  Drilled  per  Production  Unit  vs  Mining  Methods 
In  general  the  results  of  the  regressions  appear  reasonably 
stable  and  comparable.  While  development  estimates  have  wide  dis¬ 
crepancies  in  method  A,  breaking  values  are  more  consistent. 

(i)  Holes  Drilled  per  Foot  of  Stope  Development 

vs  Stope  Classification  -  Figures  51  -  55 

Except  for  Method  A  it  appears  that  about  3.5  holes 
are  drilled  and  blasted  per  foot  of  subdrift  produced. 

If  one  assumes  that  the  normal  subdrift  is  drilled  to  an 
8  foot  length  and  breaks  7.5  feet,  one  derives  a  calcu¬ 
lated  value  of  26.2  holes  per  round  and  28.0  feet  of  blast 
hole  per  foot  of  advance.  These  figures  might  appear  high 
for  the  normal  4'  x  6*  subdrift  sections  but  possibly  close 
observation  will  find  them  in  agreement  with  actual  mine 
practices . 

The  anomaly  in  holes  drilled  per  foot  of  subdrift  in 
the  "A"  method  is  understandable  when  one  considers  that 
the  labour  factors  in  this  method  were  also  quite  different 
from  other  methods.  If  much  more  stoping  is  to  be  done  with 
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the  A  mining  method  it  would  be  well  to  investigate  the  work 
and  reporting  methods  for  development  work  in  this  type  of 
st oping. 

(ii)  Holes  Drilled  per  Ton  of  Ore  Broken 

vs  Stope  Classification  -  Figures  56  -  59 

The  results  here  are  fairly  regular  and  the  A 
method  is  comparable  to  the  other  classification.  Use  factors 
in  the  narrow  widths  show  some  gaps  but  these  can  be  expected 
as  narrow  widths  are  seldom  found  in  C  type  s topes.  It 
appears  that  0.3  holes  drilled  and  blasted  per  ton  broken 
is  normal  for  ore  in  4'  -  6'  widths. 

The  6'  -  10'  width  gives  a  continuous  curve  with  an 
average  drilling  factor  of  about  0.18  holes  per  ton.  If 
one  had  a  measure  of  the  length  of  holes  drilled  it  is 
possible  some  of  the  fluctuations  would  be  removed. 

In  wide  breaking  ores  only  the  major  mining  methods  are  repre¬ 
sented.  Factors  of  0.13  and  0.10  holes  per  ton  broken 
would  be  reasonable  estimates  for  breaking  in  10'  -  20* 
and  +20'  widths. 

(2)  Pounds  of  Explosive  Reported  per  Production  Unit 

vs  Stope  Class 

In  general  the  overall  aspect  of  these  plots  is  one  of 
fluctuation.  If,  however,  one  disregards  the  B  Fay  class,  some 
measure  of  stability  is  apparent. 
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(i)  Pounds  of  Explosive  Reported  per  Foot  of 
Stope  Development  vs  Stope  Classification 

Figures  60  -  63 

Distance  from  the  mill  hole  should  have  no  effect  on 
explosive  usage.  Examination  of  the  plotted  estimates  show 
that  no  consistent  difference  is  present.  Ignoring  B  Fay 
class  we  see  that  13  lbs.  of  explosive  per  foot  of  subdrift 
is  a  reasonable  estimate  for  explosive  use  in  development. 

This  agrees  very  well  with  the  drilling  factor,  as  the  normal 
load  of  Amex  in  1-1/4"  diameter  holes  is  0.5  lbs.  per  foot. 

(ii)  Pounds  of  Explosive  Reported  per  Ton  of  Ore 
Broken  vs  Stope  Classification  -  Figures  64  -  66 

The  B  type  of  stoping  still  shows  anomalous  results. 
Ignoring  these  estimates  it  appears  that  0.8,  0.6  and  0.5  are 
reasonably  good  estimates  of  powder  use  in  widths  of  6'  -  10*, 
10'  -  20'  and  +20'.  These  values  also  are  in  good  agreement 
with  drilling  factors  when  the  normal  10'  blast  hole  is 
assumed  for  rock  breaking. 


GENERAL  DISCUSSION  OF  RESULTS 


Examination  of  the  estimates,  both  from  overall  and  from 
classified  correction  factor  regression  analyses  of  stope  contract 
data,  demonstrates  that  reasonable  and  usable  costs  and  contract  rates 
can  be  obtained  from  these  contract  records. 

Of  the  various  classifications,  it  appears  that  the  time  series 
classification  gives  the  most  useful  results.  The  mining  method 
classification  could  be  quite  useful  in  detailed  investigations  of 
particular  work  methods.  The  incentive  rate  classification  has  little 
analytical  value  except  perhaps  to  demonstrate  one  aspect  of  rate  equity. 

The  factors  calculated  with  the  regressions  of  expendable 
materials  could  be  invaluable  in  production  control.  The  use  of  this 
method  in  evaluating  the  results  of  testing  programs  on  new  materials 
and  methods  will  certainly  give  much  better  information  for  decision 
than  the  normal  raw  average  methods. 

The  regression  program  (12)  with  its  available  transformations 
has  been  a  very  useful  tool  in  this  analysis.  It  does,  however,  have 
one  major  failing.  A  number  of  the  production  units  in  the  data  were, 
under  some  regression  models,  equivalent.  All  subdrifting, for  instance, 
in  the  materials  regressions  could  have  been  considered  together  as 
one  variable  rather  than  four.  Multiple  regression  with  these 
variables  added  together  may  have,  in  this  case, returned  much  better 
results.  This  type  of  additive  transformation  capability  would  make 
this  regression  program  much  more  general  in  many  analytical  problems. 
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When  gathering  data  one  is  often  uncertain  as  to  the  best  types 
and  categories  of  measurements  to  make.  With  an  additive  transfor¬ 
mation  capability  in  a  regression  program  one  can  make  a  wide  spectrum 
of  classes  and  categories  of  measurements  with  the  assurance  that 
combination  of  like  classes  is  possible  on  analysis. 

Throughout  the  proceeding  analysis  no  mention  has  been  made 
of  possible  extensions  of  the  regression  equation  into  polynomial  or 
logarithmic  terms.  Logarithmic  terms  of  both  the  dependent  and 
independent  variables  were  tested  extensively.  Computer  outputs 
of  these  programs  are  on  file  (23,  24)  but  the  analysis  is  not 
presented  in  detail  in  this  thesis.  It  appeared  that  the  model 
Y  =£b j  Xj  +^b  1°8  Xj  improved  prediction  of  contract  hours  by  one  to 

two  percent  over  the  simple  linear  model  but  the  added  difficulty  in 
interpretation  of  the  results  far  outweighed  the  added  accuracy. 

A  second  factor  which  is  also  concerned  with  extensions  on 
the  model  is  that  of  interactions.  With  thirty-two  linear  independent 
variables  it  would  be  surprising  if  some  of  the  possible  thirty-two 
factorial  interactions  were  not  significant.  Thorough  investigation 
of  interaction  factors  could  be  part  of  a  full  systems  analysis  of 
mine  stope  production  but  would  have  only  confused  this  analysis  of  a 
relatively  simple  contract  system. 

Correlation,  a  statistical  technique  which  normally  accompanies 
the  use  of  regression,  has  been  absent  throughout  this  work.  The 
reason  for  this  exclusion  was  dictated  by  the  nature  of  the  data. 
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Correlation  is  defined  only  for  random  variables.  If  one  speaks  of 
the  statistical  correlation  of  subdrifting  and  ore  breaking,  as 
measured  from  contract  records,  then  one  implies  that  the  amount  of 
subdrifting  and  the  amount  of  ore  breaking,  done  in  a  stope  in  a  con¬ 
tract  period  is  mainly  governed  by  random  causes.  This  assumption 
immediately  negates  such  governing  factors  as  production 
and  standard  practices  which  certainly  are  operative  and  non-random 
in  the  real  situation.  Mine  contract  systems,  the  central  concern 
of  this  thesis,  are  operated  and  administered  under  real  conditions. 

For  this  reason  the  simple  correlations  which  the  computer  program 
(12)  calculates  for  all  regression,  a  sample  of  which  is  given  in 
Appendix  2,  pages  2  to  5 ,  were  not  applied  in  the  analysis. 

The  final  but  by  no  means  minor  topic  of  discussion  is  one  of 
economics.  The  following  estimate  of  the  operating  costs  of  a 
monthly  regression  analysis  of  stope  contract  records  has  been  derived 
from  the  author’s  experience  with  contract  records  and  computer 
systems  gained  in  this  work. 

The  gathering  of  the  data  in  a  form  acceptable  to  computer 
analysis  could,  for  very  little  cost,  be  incorporated  in  the  normal 
mine  accounting  system.  If  this  system  were  machine- computation 
cij-encated  the  required  data  should  be  a  byproduct  of  contract 

j 

calculations.  Computer  time  for  analysis  of  a  program  will  be  between 
5  to  10  minutes  for  a  large  machine  such  as  the  IBM  7040  and  possibly 
20  to  40  minutes  for  a  smaller  machine  such  as  the  IBM  1130.  The 
preparation  of  the  results  from  the  computer  should,  on  a  regular 
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basis,  occupy  a  clerk  no  more  than  one  to  two  hours.  The  costs  for 
the  above  should  be  no  greater  than  $20.00  per  computation  plus 
possibly  $20.00  allocated  to  clerical  workers.  If  one  compares 
this  estimate  with  the  probable  costs  of  time  study  staff  that  would 
be  required  to  produce  similar  results,  the  economic  advantages  of 
multiple  regression  are  immediately  apparent. 


CONCLUSIONS 


The  use  of  multiple  regression  as  an  analytical  and  control 
tool  for  mine  contract  administration  has  been  demonstrated.  The 
application  of  this  tool  by  a  competent  and  interested  staff  can 
provide  precise,  economic  and  timely  information  to  mine  management. 

This  information  is  the  basis  upon  which  all  good  administration  is 
founded.  One  must  not,  however,  conclude  that  this  or  other 
statistical  techniques  are  a  replacement  for  time  and  method  study 
or  the  seasoned  judgement  of  contract  administrators.  Statistical 
methods  are  tools  which,  when  applied  properly,  will  allow  the 
contract  staff  to  make  much  more  efficient  use  of  their  particular 
capabilities . 

The  techniques  developed  in  this  thesis  need  not  be  restricted 
to  underground  mining.  Contract  systems  based  on  time  study  and 
related  methods  of  evaluating  the  labour  content  of  production  are 
common  throughout  industry.  Past  records  for  payment  of  these 
contracts  should  also  be  available.  With  multiple  regression  using 
the  Correction  Factor  approach,  these  records  can  be  analysed  to 
provide  an  independent  and  unbiased  evaluation  of  most  contract  systems. 
I  feel  that  this  type  of  evaluation  would  be  an  invaluable  tool  with 
which  to  construct  better  systems,  for  all  concerned  with  contracts. 
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FIG.  41  -  CORRECTION  FACTOR  FOR  SUB  DRIFTING  50'  -  100'  FROM  MILLHOLE  VS  STOPE  CLASS 
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FIG.  42  -  CORRECTION  FACTOR  FOR  SUB  DRIFTING  100 [  -  150'  FROM  MILLHOLE  VS  STOPE  CLASS 
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FIG.  44  -  CORRECTION  FACTOR  FOR  BREAKING  +20'  WIDTHS  VS  STOPE  CLASS 


FIG.  45  -  CORRECTION  FACTOR  FOR  ELECTRIC  SLUSHING  VS  STOPE  CLASS 
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FIG.  46  -  CORRECTION  FACTOR  FOR  SECONDARY  SLUSHING  VS  STOPE  CLASS 
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COMPUTER  OUTPUT  OF  STEPWISE 
MULTIPLE  REGRESSION  ON  ALL  DATA 
UNDER  MODEL  I 


Program  No.  3 
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MEANS, VARIANCES  AND  STANDARD  DEVIATIONS 


3 

3 

3 

3 

3 

3 

3 

3 

3 

3 


1 

0. 56Q78430E 

01 

0.320237256 

03 

2 

0. 12  126666E 

02 

0. 34246257E 

0  3 

3 

0 . 90207056E 

01 

0 . 280  38  56  76 

03 

4 

0. 31380392E 

01 

0.12767206E 

03 

5 

0. 70823529E 

00 

0.2176370  f  6 

02 

6 

0.2242 9 804 E 

01 

0.922 18966E 

02 

7 

0.2568 784 3E 

02 

0.83655349E 

04 

8 

0.50752863E 

03 

0. 1 67  46036E 

06 

9 

0.2462 6 98 OE 

03 

0. 179087276 

06 

10 

0.22984863E 

0  3 

0. 3299628  IE 

06 

1 1 

0. 18232157 E 

02 

0.105295336 

05 

12 

0.8627 7 646E 

02 

0. 7227791  IE 

05 

13 

0.89831685E 

03 

0. 69683755E 

06 

14 

0.98293 3 32E 

02 

0.703  7247  9E 

05 

15 

0.21281490E 

03 

0.42443 6636 

06 

16 

0*357231376 

03 

0. 112 10 18  7 E 

07 

17 

0.3767371 26 

01 

0.616864 74E 

02 

18 

0. 179 9 6858 E 

01 

0.262  1344  76 

02 

19 

0.81021960E 

02 

0.418  9461  16 

05 

20 

0. 1226007  8G 

03 

0. 134  36107E 

06 

21 

0  •  26125490E 

01 

0.3170376  7 £ 

02 

22 

0. 1040 0  78 4 E 

02 

0.  14291336E 

04 

23 

0.26103529E 

02 

0. 1 1855544E 

04 

24 

0. 32015686E 

01 

0.861 6891 3E 

02 

25 

0. 163725496 

01 

0.682656816 

02 

26 

0. 15866666E 

02 

0. 37665379E 

04 

27 

0.24 32564 7E 

03 

0.213801726 

06 

28 

0.81741 175E 

02 

0.797400516 

05 

29 

0.  161  17647E 

00 

0.147826756 

02 

30 

0.41785723E 

01 

0. 142766876 

04 

31 

0.141 17647  E- 

"01 

0.254117646 

00 

32 

0. 46652782E 

01 

0. 56727660E 

03 

33 

0.44642352E 

00 

0.357 7 726 9E 

02 

34 

0.40420595E 

03 

0.351  3860  OE 

05 

35 

0.452392156 

01 

0.5254333 26 

01 

36 

0.24  1621  186 

03 

0. 2 1569 366 E 

05 

37 

0.713411756 

01 

0.501054996 

03 

38 

0.9636470 56 

03 

0.346850156 

06 

39 

0.20267059E 

03 

0.421 806766 

05 

40 

0.796313716 

03 

0. 288  73446E 

06 

4  1 

0. 18564377E 

01 

0.735124756 

00 

42 

0.701101696 

01 

0. 135121166 

03 

0* 1789517 46  02 
0.18505744E  02 
0 •  1674472  1C  02 
0.I1299206E  02 
0.46651588E  01 
0 . 960307  06E  01 
0.91 463298E  02 
0.4092 1921E  03 
0.42318704E  03 
0 • 5 744  2  3 89F  03 
0.102613516  03 
0.26884551E  03 
0.834767966  03 
0.26527812E  03 
0.65148801E  03 
0. 10587817E  04 
0*  78  54073  7E  01 
0.51 199070E  01 
0.20468173E  03 
0 • 36655296E  03 
0 • 56  306096E  01 
0.37803883E  02 
0. 344318816  02 
0.92827212E  01 
0 • 82  62  304  BE  01 
0.61 37212 6E  02 
0.46238699b  03 
0.28238281E  03 
0. 38448244E  01 
0. 377845036  02 
0.50410082E  00 
0.2381 7569E  02 
0.59814103E  01 
0. 1874529 3E  03 
0.2292233 2E  01 
0. 1468651 3E  03 
0.22  3842  5  7E  02 
0.58893985E  03 
0. 20537935E  03 
0.53  734017E  03 
0.857394166  00 
0. i 1 624163E  02 
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REGRESSION  COEFFICIENTS  AND  ANALYSIS  OF  VARIANCE 


3 

3 

v _ 3 

1 

1 

I 

I  13 

0 • 2648  3  3926 
0. 155147976 

03 

00 

0. 45500550 E- 02 

• 

►— 

o 

-6 

-vj  -vl 

9  9 

-4 
^  -C' 

f  3 

I 

l 

0 • 2136944  IE 

06 

0. 2 136944  IE  00 

1162.68 

0.00000 

3 

1 

1273 

0.233971356 

08 

0.183795256  05 

3 

1 

1274 

0. 44766576E 

08 

3 

1 

0. 135571  1  IE 

03 

3 

2 

0.26340110b 

03 

3 

2 

1  13 

0. 138499726 

00 

0.439808136-02 

31.49 

47.74 

3 

2 

2  1 

0.28812699E 

01 

0.205160196  00 

14.04 

7.02 

3 

2 

54.75 

3 

2 

2 

0.245103366 

08 

0.122551686  08 

769.57 

0. 00000 

3 

2 

1272 

0.202562406 

08 

0.159247176  05 

3 

2 

1274 

0.44766576E 

08 

3 

2 

0. 12619317E 

03 

0.314089506  07 

197.23 

3 

2 

0. 92735 1 1 7E 

00 

3 

3 

0.22473987E 

03 

3 

3 

1  13 

0. 122152286 

00 

0.43784  722E-02 

27.90 

47.74 

3 

3 

2  1 

0.293929386 

01 

0.194374796  00 

15.12 

7.  02 

3 

3 

3  8 

0.  104459896 

00 

0 . b61 7 5 56 BE- 02 

12.12 

4.69 

3 

3 

59.44 

3 

3 

3 

0 • 266094336 

08 

0.886981126  07 

620.89 

0.00000 

3 

3 

1271 

0. 18  15714  26 

08 

0. 142  85714E  05 

3 

3 

1274 

0.447 665 7 6E 

08 

3 

3 

0. 1  1952286E 

03 

0.209909806  07 

146.94 

3 

3 

0.83617346E 

00 

3 

4 

0.220949116 

03 

3 

4 

1  13 

0.  1153597 7E 

00 

0.423 141006-02 

27.26 

47.74 

3 

4 

2  1 

0.29428729E 

01 

0.185828756  00 

15.84 

7.02 

3 

4 

3  8 

0. 100364446 

00 

0.824709396-02 

12.17 

4.69 

3 

4 

4  20 

0.9747721 96- 

-01 

0.887639706-02 

10.98 

3.52 

3 

4 

62.96 

3 

4 

4 

0.281  84066E 

08 

0.704601666  07 

539.63 

0.00000 

3 

4 

1270 

0. 16582509E 

08 

0.130570946  05 

3 

4 

1274 

0 • 447665766 

08 

3 

4 

0. 1 1426764E 

03 

0.157463276  07 

120.60 

3 

4 

0  •  8  0  9  5  1 8  0  8  E 

00 

3 

5 

0. 200481 69E 

0  3 

3 

5 

1  13 

0. 13195886E 

00 

0.43353184E-02 

30.44 

47.74 

3 

5 

2  1 

0. 284225976 

01 

0. 178  19495E  00 

15.95 

7.02 

3 

3 

3  8 

0.8872186 46 

-01 

0.797109556-02 

11.13 

4.69 

3 

5 

4  20 

0.  10258152E 

00 

0.851 325466-02 

12.05 

3.52 

3 

5 

5  12 

0. 13226603E 

00 

0.122831816-01 

10.77 

3.10 

3 

5 

66.06 

3 

3 

5 

0.295723906 

08 

0.591447796  07 

493.97 

0.00000 

3 

5 

1269 

0. 1519418 6 E 

08 

0.119733536  05 

3 

5 

12  74 

0.44766576E 

08 

3 

5 

0. 109422826 

03 

0.138832376  07 

115.95 

3 

5 

0 . 69  76628  9E 

00 

. 
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3 

6 

0. 180768886 

03 

3 

6 

1  13 

0. 13816832E 

00 

0.424609466—02 

32.54 

47.74 

3 

6 

2  1 

0.266111056 

01 

0.173470616  00 

15.34 

7.02 

3 

6 

3  8 

0. 794942966- 

01 " 

0. 777403736-02 

10.23 

4.69 

3 

6 

4  20 

0. 1094  8958E 

00 

0. 82688719E-02 

13.24 

3.52 

3 

6 

5  12 

0. 13630333E 

00 

0.1 18909126-01 

1 1  .46 

3.10 

3 

6 

6  2 

0. 1538 1870E 

01 

0.164077976  .00 

9.37 

2.20  • 

3 

6 

68.26 

3 

6 

6 

0. 30557245E 

08 

0.509287416  07 

454.47 

0.00000 

3 

6 

1268 

0. 14209331E 

08 

0.  1  12Q6097E  05 

3 

6 

12  74 

0. 44766576E 

08 

3 

6 

0. 10585885E 

03 

0 • 984  8 5449  E  06 

67  .89 

3 

6 

0 . 66  560549E 

00 

3 

7 

0. 17635963E 

03 

3 

7 

1  13 

0.  132075446 

00 

0.4 1823 076E-02 

31  .58 

47.74 

3 

7 

2  1 

0  •  25046565E 

01 

0. 169  43438E  00 

14.78 

7.02 

3 

7 

3  8 

0  •  8420082  2E- 

-01 

0. 757002566-02 

11.12 

4.69 

3 

7 

4  20 

0. 10003805E 

00 

0.8103532 16-02 

12.34 

3.  52 

3 

7 

5  12 

0. 1374 32 926 

00 

0.  1  1550484  E-01 

i  1  .90 

3.10 

3 

7 

6  2 

0.158169276 

01 

0.159447736  00 

9.92 

2.20 

3 

7 

7  16 

0.249628456- 

“01 

0.284455826-02 

8.78 

1.82 

3 

7 

70.08 

3 

7 

7 

0.313714416 

08 

0.44816343E  07 

423.90 

0.00000 

3 

7 

1267 

0. 133951356 

08 

0.105  72325E  05 

3 

7 

1274 

0.447665766 

08 

3 

7 

0. 10282  18  IE 

03 

0.814195 9 9E  06 

77.01 

3 

7 

0.608941286 

00 

3 

8 

0. 16774834E 

03 

3 

8 

1  13 

0. 1 3004244E 

00 

0.40867855E-02 

31.82 

47.  74 

3 

8 

2  1 

0.225179056 

01 

0. 168 16272E  00 

13.39 

7.  02 

3 

8 

3  8 

0.91773612E- 

-01 

0. 7441 9127 E- 02 

12.33 

4.69 

3 

8 

4  20 

0.8094  13066- 

-01 

0. 824  6 1 32  3E-02 

9.82 

3.52 

3 

8 

5  12 

0. 14236058E 

00 

0. 1 128 1808E-01 

12.62 

3.  10 

3 

8 

6  2 

0. 1 5965994E 

01 

0.155524136  00 

10.27 

2.20 

3 

8 

7  16 

0.29  6  90715  E- 

-01 

0, 28348197E-02 

10.47 

1.82 

3 

8 

8  15 

0.379617886- 

“01 

0.46 7 55 94 4 E— 02 

8.12 

1.48 

3 

8 

• 

' 

71.56 

3 

8 

•  8 

0. 320344026 

08 

0 • 400  4300 3 6  07 

398.16 

.  0.00000 

3 

8 

1266 

0.127321736 

08 

0.10057009E  05 

3 

8 

1274 

0.4  47  665766 

08 

3 

8 

0.100284646 

03 

0.662961746  06 

65.92 

3 

8 

.  0.5 180698  IE 

00 

APPENDIX  1 


PAGE  4 


3 

9 

0. 158281 4 5E 

03 

3 

9 

1 

13 

0. 13180383E 

00 

0.400  18955  E— 02 

32.94 

47.  74 

3 

9 

2 

1 

0 • 2263289  IE 

01 

0. 16440853E  00 

13.77 

7.02 

3 

9 

3 

8 

0.80554  7706- 

01 

0. 74192197E-02 

10.86 

4.69 

3 

9 

4 

20 

0.850313048— 

01 

0.80791014E— 02 

10.52 

3.52 

3 

9 

5 

12 

0. 1458582  3E 

00 

0. 11033794E-01 

13.21 

3.  10 

3 

9 

6 

2 

0.  16509854E 

01 

0.  152  20900  E  00 

10.85 

2.20 

3 

9 

7 

16 

0. 30010040E- 

01 

0 ®  2  77 1 727 9E-02 

10.83 

1.82 

3 

9 

8 

•  15 

0. 39071050E- 

01 

0.45732835E-02 

8.54 

1.48 

3 

9 

9 

3 

0.  12970404E 

01 

0.  1 680280 1 E  00 

7.72 

1.28 

3 

9 

72.84 

3 

9 

9 

0. 32607153E 

00 

0 • 362  30 1 7  0E  07 

376.92 

0.00000 

3 

9 

1265 

0. 12  159422E 

08 

0.96121913E  04 

3 

9 

1274 

0.44766576E 

08 

3 

9 

0.9804 1784E 

02 

0.57275125E  06 

59.59 

3 

9 

0.493767Q2E 

00 

3 

10 

0. 15232350E 

03 

3 

10 

1 

13 

0. 12059339E 

00 

0 . 422  08  63  5  E~  02 

28.57 

47.  74 

3 

10 

2 

1 

0.2316 5 488E 

01 

0* 1613652 7E  00 

14.36 

7.02  „ 

3 

10 

3 

8 

0.89334326E- 

-01 

0.73755742E-02 

12.11 

4.69 

3 

10 

4 

20 

0.81 92298  3E- 

-01 

0.79329685E-02 

10.33 

3.  52 

3 

10 

5 

12 

0. 13933398E 

00 

0. 108 60 8 96 E- 01 

12.83 

3.  10 

3 

10 

6 

2 

0.  1  7040846E 

01 

0. 14941659E  00 

1 1  .40 

2.20 

3 

10 

7  16 

~0.28526147£- 

-OT' 

0.27253581E-02 

10.47 

1.82 

3 

10 

8 

15 

0 . 365430 1  IE- 

-01 

0.449 76232E-02 

8.12 

1.48 

3 

10 

9 

3 

0. 13037682E 

01 

0. 16474731E  00 

7.91 

1.28 

3 

10 

10 

9 

0. 51070582E- 

-01 

0.70873220E-0 2 

7.21 

1.07 

3 

10 

73.91 

3 

10 

10 

0. 33086952E 

08 

0 . 330 86952E  07 

358.08 

0.00000 

f 

3 

10 

1264 

0. 1 1 679624E 

08 

0 • 9240208 5 E  04 

3 

10 

1274 

0.44766576E 

08 

3 

10 

0.96 1 26003E 

02 

0.47979850E  06 

51.93 

3 

10 

0. 3981120RE 

00 

3 

1 1 

0. 15039749E 

03 

3 

1 1 

1 

13 

0. 97801385  E- 

-01 

0.511 82473E-02 

19.11 

47.74 

3 

1 1 

2 

1 

0.22  192044E 

01 

0*  158  43994E  00 

14.01 

7.02 

3 

11 

3 

8 

0. 1089238GE 

00 

0.767  10243  E-02 

14  .20 

4.69 

3 

11 

4 

20 

0.84549  147E- 

-01 

0.777 10670E-02 

10.88 

3.52 

3 

1 1 

5 

12 

0.  125061 7 8E 

00 

0. 107  95742E—01 

11.58 

3.  10 

3 

1 1 

6 

2 

0.  16730848E 

01 

0. 14627792E  00 

1 1.44 

2.20 

3 

11 

7 

16 

0.27180920E- 

-01 

0. 26730142 E— 02 

10.17 

1.82 

3 

1 1 

8 

15 

0. 37770453E- 

-01 

0.44044 1 50 E—  02 

8.58 

1.48 

3 

1 1 

9 

3 

0. 13233793E 

01 

0. 1612438 7E  00 

8.21 

1.28 

3 

1 1 

10 

9 

0.67493507E- 

-01 

0. 72695825E-02 

9.28 

1.07 

3 

1 1 

11 

10 

0.44  7  91187  E- 

-01 

0. 59396546 E— 02 

7.54 

1.  12 

3 

1 1 

75.03 

3 

1 1 

11 

0. 335901 76E 

08 

0 «  305  3652 3E  07 

345.08 

0.00000 

3 

1 1 

1263 

0. 11 176400E 

08 

0. 88490892E  04 

3 

11 

1274 

0 o  44766576E 

08 

3 

1 1 

0.94069598E 

02 

0. 50322375E  06 

56.87 

3 

11 

0 • 23754  56  7E 

00 
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3 

12 

• 

0.  1427696 1 E  03 

3 

12 

1 

13 

0.91561695 E— 01 

0.501995066-02 

18.24 

47.74 

3 

12 

2 

1 

0.21585216E  01 

0.1  539967 3E  00 

'  14.02 

7.02 

3 

12 

3 

8 

0. L  18291986  00 

0. *'7 52  3963 6 E- 02 

15.72 

4.  69 

3 

12 

4 

20 

0. 774037436-01 

0.758904546-02 

10.20 

3.52 

3 

12 

5 

12 

0. 11885  7206  00 

0. 10506122E-01 

11.31 

3.10 

3 

12 

6 

2 

0.161826376  01 

0.142 1697 3E  00 

11.38 

2.20 

3 

12 

7 

16 

0.285624976-01 

0. 2 60 01868 E- 02 

10.98 

1.82 

3 

12 

8 

15 

0.367634776-01 

0. 42781426 E- 02 

8.59 

1.48 

3 

12 

9 

3 

0.130950806  01 

0. 15657310E  00 

8.36 

1.28 

3 

12 

10 

9 

0. 76993184 E- 01 
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32 

18 

11 

0*  13500177E  00 

0. 21606362 E- 01 

6.25 

0.79 

3 

32 

19 

28 

0. 67872708 E— 01 

0.88A010A6E-02 

7.68 

0.65 

3 

32 

20 

6 

0.  1A973139E  01 

0.22930307E  00 

6.53 

0.65 

3 

32 

21 

27 

0.25075913 E— 01 

0.50100027E-02 

5.01 

0.  3 A 

3 

32 

22 

5 

0.2A59905AE  01 

0.A9AA0731E  00 

A. 98 

0.29 

3 

32 

23 

30 

0*221  50695E  00 

0.5886911AE-01 

3.76 

0.  18 

3 

32 

2  A 

26 

0.  12986A20E  00 

0. 358 85877 E— 01 

3.62 

0.17 

3 

32 

25 

2A 

0.7A208A71E  00 

0. 235 A5036E  00 

3.15 

0.  1A 

3 

32 

26 

22 

0.23A76116E  00 

0*  7666601 A E— 01 

3.06 

0.09 

3 

32 

27 

21 

0 • 92 A 5A668E  00 

0 . 3869098 2 E  00 

2.39 

0.09 

3 

32 

28 

25 

0*  550527A2E  00 

0 • 2789309 1 E  00 

1.97 

0.07 

3 

32 

29 

23 

0  •  169AA7A8E  00 

0.823411A0E— 01 

2.06 

0.  05 

3 

32 

30 

32 

0.160A9157E  00 

0. 917 10 66 8 E— 01 

1.75 

O.OA 

3 

32 

31 

29 

0.8272A985E  00 

0*552  53606E  00 

1.50 

0.03 

3 

32 

32 

33 

— 0.9973288AE— 01 

0.35998A35E  00 

-0.28 

0.00 

3 

32 

8A.  A5 

3 

32 

32 

0.3780A76AE  08 

0*118 13989E  07 

210.76 

0.00000 

3 

32 

12A2 

0*696 1811AE  07 

0.56053232E  0A 

3 

32 

127A 

0 • AA 766576E  08 

3 

32 

0*  7A868706E  02 

0* A3000000E  03 

0.08 

3 

32 

0.6A3A8152E-0 2 
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MULTIPLE  REGRESSION  ON  ALL  DATA 
UNDER  MODEL  II 


Program  No.  60 
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MEANS, VARIANCES  ANO  STANDARD  DEVIATIONS 


60 

1 

0.807  302906 

01 

0. 77932906E 

03 

0. 2791 6466E 

02 

60 

2 

0.83673833E 

02 

0. 16304630E 

05 

0.  12768958E 

03 

60 

3 

0 .66  7  53 1 2  5E 

02 

0. 15353902E 

05 

0. 12391086E 

03 

60 

4 

0 • 24  790499E 

02 

0. 7968009  IE 

04 

0. 89263705E 

02 

60 

5 

0 • 59  49 1 75  7E 

01 

0. 15356470E 

04 

0.39  187332E 

02 

60 

6 

0. 14355068E 

02 

0.37772871E 

04 

0.61 45963  7E 

02 

60 

7 

0. 17210847E 

02 

0 • 375  52875E 

04 

0.612804006 

02 

60 

8 

0 . 24  36 1 30  3E 

03 

0. 38583035E 

05 

0. 196425656 

03 

60 

9 

0. 10343319E 

03 

0. 31590975E 

05 

0.  17773850E 

03 

60 

10 

0  •  89640904E 

02 

0.501 87 32 2E 

05 

0.22402527E 

03 

60 

11 

0. 10939293E 

02 

0. 37906316E 

04 

0 • O 1 568 1 05E 

02 

60 

12 

0. 4313882 3E 

02 

0.  1 80  6947  8  E 

05 

0. 13442276E 

03 

60 

13 

0. 35932567E 

03 

0. 11149431E 

06 

0. 33390764t 

03 

60 

14 

0. 19658650E 

02 

0.28148976E 

04 

0.5  3  05  560  8E 

02 

60 

15 

0 . 25537776E 

02 

0.611 18853E 

04 

0. 78  178547E 

02 

60 

16 

0. 357231 1  IE 

02 

0. 112 10184E 

05 

0. 10587816E 

03 

60 

17 

0.565  10563E 

02 

0.  1  38  79456E 

05 

0.  1 1  78 1 1 1  Oh 

03 

60 

18 

0.24295760E 

02 

0 • 477  7400  8  E 

04 

0.69118743E 

02 

60 

19 

0 • 24306578E 

02 

0. 37705137E 

04 

0.61404508E 

02 

60 

20 

0 . 3678022  7E 

02 

0.  1 20  9249  5E 

05 

0. 10996588E 

03 

60 

21 

0. 195941 18E 

01 

0. 17833369E 

02 

0.42229573E 

01 

60 

22 

0.20801568E 

02 

0.571 65343E 

04 

0. 75607765E 

02 

60 

23 

0 . 326294 1 2E 

02 

0.  185  2428  8  E 

04 

0. 4303985  IE 

02 

60 

24 

0 . 800392 1 5E 

01 

0. 53855570E 

03 

0. 2320680  3E 

02 

60 

25 

0.49117647E 

01 

0.61439114E 

03 

0.24786915E 

02 

60 

26 

0 • 6346662  3E 

01 

0 • 602  6455  8E 

03 

0. 2454884  IE 

02 

60 

27 

0 • 364884  3  1  E 

02 

0.481 0538  3E 

04 

0.69358044E 

02 

60 

20 

0. 31061633E 

02 

0. 115 14459E 

05 

0. 10730545E 

03 

60 

29 

0. 161  17647E 

00 

0. 147  8267  5E 

02 

0.38  448244E 

01 

60 

30 

0.41705723E 

01 

0. 14276687E 

04 

0. 37784503E 

02 

60 

31 

0. 141176476-01 

0.2541 1764E 

00 

0. 50410082E 

00 

60 

32 

0.46652782E 

01 

0.567276606 

03 

0. 23817569E 

02 

60 

33 

0  •  44642352E 

00 

0. 35777269E 

02 

0.598141036 

01 

60 

34 

0. 14561471E 

04 

0.456034216 

06 

0.675303046 

03 

60 

35 

0.452 39215E 

01 

0 • 52543332E 

01 

0.229223326 

01 

60 

36 

0.241621 18E 

03 

0.215693666 

05 

0. 1468651 3E 

03 

60 

37 

0.71341 175E 

01 

0.501 05499E 

03 

0.223842576 

02 

60 

38 

0.96364705E 

03 

0.346850156 

06 

0.588939856 

03 

60 

39 

0.20267059E 

03 

0.421806766 

05 

0. 20537935E 

03 

60 

40 

0.79631371E 

03 

0.288  73446E 

06 

0. 53  73401 7  E 

03 

60 

41 

0.185643776 

01 

0.735 12475E 

00 

0.857394166 

00 

60 

42 

0.701101696 

01 

0.135121166 

03 

0.116241636 

02 
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SIMPLE  CORRELATION  COEFFICIENTS 


60 

I 

2 

0.07094 

-0.02720 

-0.02048 

-0.01324 

0.01011 

0.03736 

60 

I 

8 

0.06176 

0.08251 

0.22  180 

-0.02877 

-0.04843 

0.26953 

60 

1 

14 

-0.02657 

0.17608 

0. 15745 

0. 12293 

0.00090 

0.16999 

\  60 

1 

20 

0.04388 

0.06711 

0.01818 

0.11247 

0.05895 

0.24497 

(  60 

1 

26 

0.08079 

0. 15788 

-0.04004 

0.01706 

-0.01526 

-0.00891 

60 

I 

32 

0.03780 

0.02785 

0.44130 

0.31485 

0.18610 

-0.00851 

60 

1 

38 

0.18116 

0.14651 

0.  13875 

-0. 00585 

0.09182 

60 

2 

3 

“0.01624 

-0.11478 

-0.06823 

0.09580 

0.03237 

0.07379 

60 

2 

9 

-0.11919 

-0.04088 

0.05459 

0.02953 

-0.11221 

-0.02472 

60 

2 

15 

-0.03183 

-0.06319 

-0.07505 

-0.04073 

-0.09205 

-0.09865 

60 

2 

21 

0.00408 

-0.09530 

-0.07246 

-0.04317 

-0.00414 

-0.03073 

60 

2 

27 

-0.00267 

0.09670 

-0.00737 

0.06453 

0.05433 

-0.00612 

60 

2 

33 

0.00517 

0.09829 

0.03457 

0.3753b 

0.15681 

0.29818 

60 

2 

39 

0.16021 

0.23495 

-0.09344 

-0.00563 

60 

3 

4 

0.08930 

-0.05163 

-0.04430 

-0.01003 

0.18130 

-0.04602 

60 

3 

10 

-0.05957 

-0.01372 

-0.00681 

-0.00295 

0.00368 

-0.06693 

60 

3 

16 

-0.03243 

-0.07307 

-0.07162 

-0.08200 

-0.06303 

0.04886 

60 

3 

22 

-0.04422 

0.02665 

-0.06005 

0.02973 

-0.04309 

-0.01038 

60 

3 

28 

-0.05059 

-0.02260 

-0.01058 

-0.01510 

-0.02559 

0.01 6^3 

60 

3 

34 

0.11100 

0.03843 

0.40060 

0.  14271 

0.34664 

0.  16843 

60 

3 

40 

0.35496 

0.03830 

0.03038 

60 

4 

5 

0.06739 

-0.01931 

-0.03048 

0. 11747 

-0.08118 

-0.03323 

60 

4 

1  1 

0.03124 

0.01406 

-0.04753 

0.  12854 

-0.07541 

0. 00863 

60 

4 

17 

-0.08875 

-0.04318 

-0.03543 

-0.02672 

0.08722 

-0.04325 

60 

4 

23 

0.00235 

-0.05400 

-0.03496 

-0.01981 

-0.00144 

-0.00127 

60 

4 

29 

-0.01165 

-0.02757 

-0.00778 

-0.01445 

-0.00489 

0.06130 

60 

4 

35 

0.03506 

0.24344 

0.09861 

0. 19918 

0.0992  t 

0. 18953 

60 

4 

41 

-0.00160 

0.02043 

60 

5 

6 

0.01406 

0.07760 

0.12638 

-0.04863 

-0.05172 

0.02086 

60 

5 

12 

0.06983 

0.00023 

0.32064 

-0.01840 

-0.05126 

-0.05338 

60 

5 

10 

-0.00097 

-0.05230 

-0.03747 

0.03012 

-0.03697 

0.02683 

60 

5 

24 

-0.00527 

-0.01878 

-0.03  152 

0.08465 

-0.00991 

-0.00637 

60 

5 

30 

-0.01242 

-0.00425 

0.01 131 

0. 18627 

0.08159 

0.05262 

60 

5 

36 

0.  14837 

0.03884 

0.15047 

0. 10264 

0.09059 

-0.00560 

60 

6 

42 

0.08797 

60 

6 

7 

0.20368 

-0.06423 

-0.02183 

-0.05269 

0.03151 

0.01482 

60 

6 

13 

-0.06105 

0.03974 

0.02819 

-0.01866 

-0.03060 

0.02182 

60 

6 

19 

0.00859 

0.02457 

-0.03432 

“0.00235 

-0.06857 

0.03862 

60 

6 

25 

-0.01288 

0.02455 

-0.00369 

0.03895 

0.06415 

0.20397 

60 

6 

31 

-0.00655 

-0.02687 

0.01797 

0.07247 

0.04149 

0.15654 

60 

6 

37 

0.07146 

0.11797 

0.11213 

0.05947 

-0.06265 

-0.00787 

60 

7 

8 

-0.05768 

-0.08333 

-0.08604 

0.02011 

0.01611 

-0.00418 

60 

7 

14 

0.00814 

0.07347 

-0.06160 

0.00892 

0.07387 

0.01404 

60 

7 

20 

0.09610 

-0.04912 

0.16296 

“0.06130 

0.16713 

-0.03779 

60 

7 

26 

0.02894 

-0.05593 

0.10172 

0.02933 

0.23705 

-0.00787 

60 

7 

32 

0.08798 

0.04000 

0. 10753 

0.08513 

0.10941 

0.03074 

60 

7 

38 

0. 12430 

0.08922 

0.05452 

0.03157 

-0.02798 

60 

8 

9 

-0.03680 

-0.06665 

-0.02184 

0.01019 

0.31269 

0.  14190 

60 

8 

15 

-0.06746 

0.00412 

-0.08876 

-0.00432 

-0.07229 

0.09642  t 

60 

8 

21 

0.11089 

0.12677 

0.24390 

-0.05931 

0.04383 

-0.04832 

60 

8 

27 

0.12320 

0.00044 

-0.01281 

-0.07512 

-0.02675 

-0.05668 

60 

8 

33 

0.02706 

0.41547 

0.15715 

0.52435 

0.06220 

0.57561 

06 
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60 

8 

39 

0. 16486 

0.57370 

0.24946 

0. 12248 

60 

9 

10 

0.01236 

-0.06363 

-0.08840 

0.37981 

-0.02345 

0.12358 

60 

9 

16 

0.  15970 

0.13460 

0.04690 

0.07935 

0.13811 

0.07814 

60 

9 

22 

0.13720 

0.18601 

0.08874 

0.09084 

0.05412 

0.13278 

60 

9 

28 

-0.07556 

-0.01407 

-0.04906 

-0.00373 

-0.00654 

-0.01232 

60 

9 

'34  0.35081 

0.25290 

0.18451 

-0.06636 

0.26332 

0. 17402 

60 

9 

40 

0.22157 

0.26439 

0 •  12  309 

60 

10 

1 1 

-0.05314 

-0.08498 

0.50708 

-0.09746 

0.01316 

0.17181 

60 

10 

17 

0.06144 

-0.08617 

0.11742 

0.01001 

0.08489 

-0.00632 

60 

10 

23 

0.52721 

0.07876 

0.18288 

-0.04317 

0.09915 

-0.06959 

60 

10 

29 

-0.01679 

-0.02946 

-0.01121 

-0.02500 

-0.00444 

0.39185 

60 

10 

35 

0.28406 

0.22410 

-0.04305 

0.26505 

0.25917 

0. 17970 

60 

10 

41 

0.34214 

0.23427 

60 

11 

12 

-0.05629 

-0.18613 

-0.04518 

-0.05308 

-0.06000 

-0.08529 

60 

11 

18 

-0.05173 

-0.06461 

-0.05657 

-0.01112 

-0.04645 

-0.07077 

60 

11 

24 

0.00890 

-0.03153 

-0.00650 

0.04148 

0.02752 

-0.00745 

\  - 

60 

11 

30 

0.03001 

-0.00498 

0.09776 

-0.01262 

-0.05500 

-0.01576 

k 

60 

11 

36 

0.04240 

0.01137 

0.00563 

0.01663 

0.00625 

-0.05311 

( 

60 

1 1 

42 

0.03221 

60 

12 

13 

-0.34322 

0.08946 

-0.09741 

-0.09805 

-0.14249 

-0.07952 

60 

12 

19 

-0.10494 

-0.10484 

-0.01 545 

-0.08806 

-0.00266 

-0.06112 

I 

60 

12 

25 

-0.05815 

-0.07510 

-0.00140 

0.49223 

-0.01346 

0.  1 045*7 

60 

12 

31 

-0.00899 

-0.02004 

0.08136 

-0.04419 

-0.08101 

-0.01058 

60 

12 

37 

-0.00630 

0.00344 

-0.03811 

0.00455 

0.05304 

0.00607 

. 

60 

13 

14 

-0.01606 

0.10113 

0.22823 

0.20864 

0.03380 

0.21492 

60 

13 

20 

0. 17288 

0.08473 

0.17646 

0.49649 

0.13233 

0.17647 

60 

13 

26 

0.02331 

0.29454 

-0.22948 

-0.00866 

-0.05714 

-0.01985 

60 

13 

32 

-0.01745 

-0.01688 

0.69091 

0.46491 

0.33625 

-0.05094 

60 

13 

38 

0.40813 

0.29658 

0.34845 

0.47023 

0.23820 

60 

14 

15 

-0.07179 

-0.05823 

-0. 10545 

-0.04700 

-0.05646 

-0.04836 

60 

14 

2  1 

0.09005 

-0.05721 

-0.00339 

-0.00771 

-0.05755 

-0.02124 

60 

14 

27 

0.13133 

0.00509 

-0.01554 

-0.01078 

-0.01039 

-0.00771 

60 

14 

33 

0.00455 

0.08669 

-0.02125 

0.  12805 

0.01962 

0. 11514 

60 

14 

39 

0.05268 

0.08763 

-0.01114 

-0.02070 

60 

15 

16 

-0. 10546 

0.21279 

0.22679 

0. 20281 

0.26393 

-0.05093 

60 

15 

22 

0.10354 

-0.08953 

0.05277 

0.09423 

0.10281 

-0.01226 

60 

15 

28 

-0.09200 

0.04512 

-0.00001 

-0.00916 

0.03062 

0.01062 

60 

15 

34 

0.21402 

0.17879 

-0.03688 

0.01848 

-0.01477 

0.02632 

60 

15 

40 

-0.00990 

0.07977 

-0.01 116 

60 

16 

17 

0.42629 

0.00519 

0.31032 

0.  16618 

0.00555 

0.05679 

i 

60 

16 

23 

0.12428 

0.01941 

0.24881 

0.01709 

0.07623 

-0.09774 

j 

60 

16 

29 

-0.00212 

0.00563 

-0.00919 

0.00035 

-0.00174 

0.31695 

r 

60 

16 

35 

0.31034 

0.04324 

-0.00591 

0.08681 

0.05690 

0. 04269 

60 

16 

41 

0.22729 

0.13142 

60 

17 

18 

0.01823 

0.33684 

0.35011 

-0.07903 

0.28497 

-0.01432 

60 

17 

24 

0.03773 

0.05979 

0.17488 

0.04738 

-0. 1322  5 

-0.02012 

60 

17 

30 

-0.00787 

-0.01344 

0.06886 

-0.02186 

0.29655 

0.29538 

60 

17 

36 

-0.09785 

0.00409 

-0.04234 

-0.00568 

-0.06217 

0.21451 

60 

17 

42 

0.04747 

60 

18 

19 

0.02222 

0.22837 

-0.05001 

0. 17835 

-0.13301 

0.05028 

60 

18 

25 

-0.03532 

0.09930 

-0.01471 

-0.09147 

0.05636 

0.00282 

60 

18 

31 

-0.00985 

0.02631 

-0.00913 

0. 12287 

0.08865 

-0. 05876 

60 

18 

37 

0.01 1‘93 

-0.01846 

-0.00712 

-0.01113 

0.02373 

0.00483 

60 

19 

20 

0.18689 

-0.00969 

0.09303 

0.02369 

0.12022 

0. 19094 
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60 

19 

26 

0. 14051 

0.01004 

-0.09515 

0.01183 

0.01056 

-0.01109 

60 

19 

32 

0. 10891 

-0.00736 

0.29086 

0.25335 

-0.07256 

-0.00817 

60 

19 

38 

-0.05603 

0.00565 

-0.07100 

0.  13043 

0.07808 

60 

20 

21 

-0.01575 

0.64044 

0.01615 

0.04853 

0.05560 

0.13443 

60 

20 

27 

0.07327 

-0.09655 

0.07475 

-0.03257 

-0.00937 

0.02450 

60 

20 

33 

-0.00313 

0.31287 

0.22875 

0.00822 

0.01933 

0.10253 

60 

20 

39 

-0.01259 

0.10299 

0.28880 

0.05805 

60 

21 

22 

-0.00806 

0.10423 

-0.01609 

0.05093 

0.01096 

0.00487 

60 
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REGRESSION  COEFFICIENTS  AND  ANALYSIS  OF  VARIANCE 
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LIST  OF  WORKING  PLACES  -  (Continued) 


Working  No.  of  Cent.  Rec. 

Place  Jan./ 63  Feb./ 62  Mining  Ore  Ground 

No.  June . / 64  Nov./64  Method  Control  Control  Breaking 
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summary  statistics  on  rf.s i duals' 


SUM  OF  RESIDUALS 

SUM  OF  SQUARES  OF  RESIDUALS 

MEAN  OF  RESIDUALS 

STANDARD  DEVIATION  OF  RESIDUALS 

DURR  IN-WATSON  STATISTIC 
CORRELATION  OF.RESIDUAL5 
NUMRFR  OF  RUNS 

PUSTTTVT~5TGN'§ 

negative  signs 

Mg  AN 

STANDARD  DEVIATION 
coefficient  OF  SKEWNESS 
STANDARD  DEVIATION 
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APPENDIX  3 


CONTROLS  ON  THE 
GENERATION  AND  CLASSIFICATION 


OF  THE  DATA 


1 


APPENDIX  3 

ELDORADO  MINING  AND  REFINING 
STOPE  CONTRACT  RATES 
February  1962  -  November  1964 


GENERAL 


All  the  following  quantities  will  be  based  on  engineering 
figures.  Tonnage  broken  will  be  calculated  on  the  basis  of  12  cu. 
ft. /ton. 


Contractors  are  required  to  maintain  timber  work  in  and  around 
their  working  places. 

Contractors  are  responsible  for  the  cleaning  up  of  backfill 
spill  when  caused  by  poor  timbering  practices. 


BREAKING 


Width  (Ft.) 

Rate/Ton 

Amex 

Allowance 
Lbs . /Ton 

High 

Explosive 

Allowance 

0-4 

1.20 

6.0 

4-6 

.67 

3.0 

6-10 

.48 

1.5 

0.33 

10  -  20 

.42 

1.1 

lbs . / fuse 

Over  20 

.39 

0.8 

Slashing 

.48 

1.5 

(except  breast  equivalent) 

Subdrifting  (Basic) 

6.90/ft. 

12.0/ft, 

3.0 

Stope  Crew  Raising 

. 

lbs ./ft 

815 

Block  Holing  &  Chute  Blasting 

lbs ./ft. 

As  required 
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BREAKING  -  (continued) 

1.  Fuse,  caps  and  Thermalite  to  be  charged  at  the 
usual  rates. 

2.  No  charge  for  powder  used  up  to  the  maximum 
allowance. 

3.  When  the  maximum  allowance  has  been  exceeded, 

the  excess  over  the  maximum  allowance  will  be 

charged  to  the  contract  at  the  following  rates: 

Amex  -  0.15/lb. 

Digel  or 

Gilgel  -  0.30/lb. 

4.  The  price  for  subdrifting  will  increase  $0. 50/ft. 
for  each  50  feet  increase  in  scraping  distance. 


SLUSHING 


Slushing 

5  H.P.  Air  Slusher 
10  H.P.  Air  Slusher 
Electric  Slusher 


Primary 


Secondary 


One  Slusher  Two  S lusher s 


$0.60  Per  Ton  $0.20  Per  Ton  $0.15  Per  Ton 

0.50  "  "  0.20  "  "  0.15  "  " 

0.40  Per  Ton  0.20  Per  Ton  0.15  Per  Ton 


STOPE  TIMBERING 

Manway  and  Millholes  (complete  with  Bulkhead) 

Stulled  (30%  each  side  and  10%  for  lining)  $13.50  Per  Ft. 
Cribbed  (30%  each  side  and  10%  for  lining)  15.00  Per  Ft. 
In  Square  Set  Stopes  5.00  Per  Ft. 


.3,5X02.0$  s* 


.  '  - 
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STOPE  TIMBERING  -  (continued) 


Gob  Fence 

$0.30 

Per  Sq.  Ft 

Hangingwall  Fill  Fence 

0.30 

Per  Sq.  Ft 

Sill  Fence 

0.25 

Per  Sq.  Ft 

Stope  Fence  or  Drift  Wall  Fence 

0.20 

Per  Sq.  Ft 

Pillar  Fence 

0.40 

Per  Sq.  Ft 

Slusher  Floor  -  0  -  7'  wide, 

14"  spacing 

0.12 

Per  Sq.  Ft 

Slusher  Floor  -  7'  up  wide, 

14"  spacing 

0.10 

Per  Sq.  Ft 

Slusher  Floor  -  less  than  14"  spacing 

0.12 

Per  Sq.  Ft 

Sill  or  Lagging  Floor 

0.18 

Per  Sq.  Ft 

Test  Holes 

0.75 

Each 

Permanent  Posts  with  Squeeze  Block 

2.50 

Each 

Crib  Back  Support  41  x  4’ 

3.00 

Per  Ft. 

Permanent  Bulkheads 

0.40 

Per  Sq.Ft. 

Bearer  Sets  for  Manway  and  Millhole 

8.00 

Each 

"A"  Sets 

8.00 

Per  Set 

Boom 

3.00 

Each 

Square  Sets 

10.00 

Each 

Bracing  Sets 

1.00 

Each 

Rockbolts  -  Hand  installed  w/Plank 

3.50 

Each  12.2 

Rockbolts  -  Hand  installed  no  Plank 

2.50 

Each 

Rockbolts  -  Machine  installed  w/Plank 

2.00 

Each 

Rockbolts  -  Machine  installed  no  Plank 

1.25 

Each 

59 
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TABLE  1 

VARIABLE  LISTING  FOR  STOPE  CONTRACT  RECORDS  DATA 


Variable 

No. 

Description 

Units  of 
Measurement 

1. 

Credit  Time 

hours 

2. 

Subdrifting  O'  -  50'  from  millhole 

feet 

3. 

Subdrifting  50*  -  100*  from  millhole 

*feet 

4. 

Subdrifting  100'  -  150'  from  millhole 

feet 

5. 

Subdrifting  +  150'  from  millhole 

feet 

6. 

Stope  raising 

feet 

7. 

Ore  broken  4'  -  6'  widths 

*tons 

8. 

Ore  broken  6'  -  10*  widths 

tons 

9. 

Ore  broken  10 '  -  20 '  widths 

tons 

10. 

Ore  broken  +  20'  widths 

tons 

11. 

Slushing  with  5  HP  air  slushers 

tons 

12. 

Slushing  with  10  HP  air  slushers 

tons 

13. 

Slushing  with  Electric  slushers 

tons 

14. 

Secondary  slushing 

*tons 

15. 

Slusher  floor  0*  -  71  widths 

^square  feet 

16. 

Slusher  floor  +  7'  widths 

^square  feet 

17. 

Cribbed  millhole 

lineal  feet 

18. 

Stulled  millhole 

lineal  feet 

19. 

Gob  fence 

*square  feet 

20. 

Hangingwall  fence 

square  feet 

21. 

Test  holes 

each 

22. 

Rockbolts  with  plates 

each 

23. 

Rockbolts  plain 

each 

24. 

Posts 

each 

25. 

Cribbing 

feet 

26. 

Bulkhead 

square  feet 

27. 

Tramming 

tons 

28. 

Muck  and  tram 

*tons 

29. 

Miscellaneous  Development 

dollars 

30. 

Miscellaneous  Ore  Breaking 

^dollars 

31. 

Miscellaneous  Slushing 

dollars 

32. 

Miscellaneous  Timber  Work 

*dollars 

33. 

Miscellaneous  Credit 

dollars 

34. 

Contract  Time 

hours 

35. 

Number  of  Workmen  on  Contract 

each 

36. 

Number  of  Blastholes  Drilled 

each 

37. 

Number  of  Ream  Holes  Drilled 

each 

38. 

Explosives  Use  Reported  by  Workmen 

pounds 

39. 

High  Explosives  as  Calculated  by  Supervision 

pounds 

40. 

Amex  as  Calculated  by  Supervision 

pounds 

41. 

Incentive  Rate 

dollars/hr . 

42. 

Penalty  Charges 

*dollars 

NOTE:  For  variables  marked  *  see  notes  taken  during 
data  preparation  Appendix  3,  page  5. 


.  ■ 

: '  ■  - 

- 


, 

■ 


NOTES  ON  CONTRACT  INFORMATION 


Major  part  of  Miscellaneous  Breaking  is  Breasting  @  $1. 20/Ton 
(-45  wide)  and  $1. 30/Ton. 

Mucking  and  tramming  @  .38$  includes  a  proportion  @  ,43c  and 
some  @  ,25c . 

Slusher  Floor  @  .12$  and  .  10c  -  price  and  small  method  change 
in  early  1963  put  most  of  floor  @  .10$.  Prior  to  this  it 
was  .12$. 

The  major  penalty  charges  less  than  $12.00  are  for  lost  bits, 
pipe  wrenches,  tools  etc.  Charges  in  excess  of  $12.00  usually 
include  some  charge  for  excess  powder  usage  -  High  explosive. 

Breaking  @  .48c  includes  all  slashing  as  well  as  6'  -  10* 
widths  breasted. 

Double  slush  @  .20c  includes  some  double  slush  @  .15c* 
Miscellaneous  timber  includes  Bearer  Timbers  @  $8.00  each. 


$1. 00/Ft. 
0.50/Ft. 

0. 15/Sq.Ft 
8.00  Each 
8.00  Each 


@ 

@ 

@ 

@ 

@ 


Manway 

Slide 


Planking 


A  sets 
Stulls 


Stopes  such  as  744  where  square  setting  is  done  are  not  in¬ 
cluded  in  the  data.  Also  much  of  964,  where  the  major  work 
has  been  for  development  for  long  hole  s toping. 


Gob  fence  @  .30c  includes  some  Sill  Fence  @  .25c/Sq.Ft. 
Some  Rse  @  $7.65/Ft.  included  in  S.  Dr.  @  $7. 40/Ft. 
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TABLE  2 

CLASSIFICATION  SYSTEMS  FOR  S TOPING  CONTRACT  RECORDS 


I.  Classification  by  Mining  Method  -  (based  on  the  judgment  of 

the  Chief  Engineer) 

A.  Standard  Hangingwall  cut  and  fill,  average  longitudinal 
10 '  cuts  at  right  angles  to  the  hangingwall  with  cribbed 
mill  holes,  hangingwall  rockbolts.  Hangingwall  fences 
and  floor  on  each  lift. 

B.  Fay  &  Verna  type  footwall  cut  and  fill;  Longitudinal  101 
cuts  at  right  angles  to  the  hangingwall.  Stulled  and 
cribbed  mill  holes,  some  hangingwall  fences,  floor  on 
each  lift.  Irregular  footwall  requires  some  subdrifting 
and  slashing. 

C.  Verna  type  cut  and  fill.  Irregular,  generally  flatter, 
dipping  orebodies  mined,  by  open  pillar  stoping  methods 
through  footwall  boxholes.  Filling  done  here  mainly  for 
ground  support  after  a  section  has  been  mined  up  to  a 
higher  boxhole.  Flooring  is  relatively  small  and  much 
gob  and  hangingwall  fences  are  used  to  close  off  filled 
areas . 

D.  Fay  West  open  stoping.  Ore  bodies  similar  to  C  but  are 
smaller,  more  erratic  and  steeper.  Mining  method  is 
similar  except  that  often  a  complete  orebody  is  extracted 
before  filling. 

E.  Open  track-level  stoping.  Used  for  isolated  ore  bodies 
and  remnants.  Ore  is  mined  by  open  stoping  with  posts 
and  pillars.  Ore  is  slushed  down  and  mucked  into  cars 
at  track  level. 

F.  Static  shrink.  Stope  is  mined  similar  to  Fay  West 
open  stopes  but  sufficient  ore  is  kept  in  the  stope 
to  provide  a  working  floor.  On  completion  of  break¬ 
ing,  the  ore  is  slushed  out  to  level  chutes  with 
rockbolts  and  posts  used  as  ground  control.  The 
empty  stope  may  or  may  not  be  filled. 
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TABLE  2  (Continued) 


II,  Classification  by  Ore  Control.  (based  on  judgment  of  Mine  Staff) 

Regular  =  1 
Good  =  2 
Erratic  =  3 


III.  Classification  by  Ground  Control,  (based  on  judgment  of  Mine  Staff) 

Normal  =  1 

Heavy 

Hangingwall  =  2 

Heavy  back  =  3 

IV.  Classification  by  Ore  Breaking  Characteristics. 

(based  on  judgment  of  Mine  Staff) 

Easy  =  1 

Normal  -  2 

Hard  =  3 


V.  Classification  by  Time  Periods. 

(Given  in  C.C.  7,  8,  9  of  all  data  records) 

Data  can  be  sorted  into  month  and  year  time  intervals. 


VI.  Classification  by  Level  of  Incentive  or  Bonus  Payments. 
(Given  in  C.  C.  61  -  70  of  data  record  6) 

Data  was  sorted  into  classes  of  $0. 50/hr.  increments. 
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LIST  OF  WORKING  PLACES 
SAMPLED  FOR  REGRESSION  ANALYSIS 
FEBRUARY  1962  -  NOVEMBER  1964 


Working 

Place 

No. 

No.  of 

Jan. / 63 
Jun. / 64 

Cont.  Rec . 

Feb . / 62 

Nov . / 64 

Mining 

Method 

Ore 

Control 

Ground 

Control 

Breaking 

0243 

12 

16 

E 

3 

1 

2 

0255 

3 

3 

E 

3 

1 

2 

0248 

1 

D 

3 

1 

1 

0253 

1 

f(s) 

3 

? 

? 

0293 

1 

e(S) 

3 

3 

2 

0304 

1 

b(S) 

1 

2 

2 

0343 

14 

31 

E 

3 

1 

2 

0355 

12 

17 

F 

3 

1 

2 

0374 

3 

C 

3 

1 

2 

0381 

2 

E 

3 

1 

2 

0385 

8 

C 

3 

1 

2 

0443 

10 

1 

3 

2 

0455 

30 

38 

D 

3 

1 

2 

0477 

15 

23 

C 

3 

1 

2 

0457 

1 

F 

1 

1 

2 

0473 

9 

E 

3 

1 

2 

0501 

3 

A 

1 

2 

2 

0543 

5 

8 

D 

1 

3 

2 

0578 

33 

43 

B 

1 

3 

2 

0581 

26 

40 

B 

1 

3 

2 

0586 

19 

32 

C 

3 

1 

2 

0593 

2 

2 

?<s) 

0597 

10 

18 

B 

? 

0594 

6 

C 

3 

1 

2 

0632 

1 

?(S) 

0606 

1 

B(S) 

1 

2 

2 

0616 

1 

B(S) 

1 

1 

2 

0672 

7 

14 

C(S) 

1 

2 

0673 

8 

13 

B 

1 

2 

0644 

1 

B 

3 

2 
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LIST  OF  WORKING  PLACES  -  (Continued) 


Working 

Place 

No . 

No,  of 

Jan. / 63 
June . / 64 

Cont.  Rec. 

Feb./62 

Nov . / 64 

Mining  Ore  Ground 

Method  Control  Control 

Breaking 

0678 

4 

B 

3 

2 

0681 

8 

B 

3 

2 

0686 

4 

B 

3 

2 

0689 

2 

B 

1 

2 

0697 

2 

B 

1 

2 

0773 

19 

67 

B 

1 

2 

0774 

2 

7 

?  (S) 

0776 

12 

17 

E 

1 

2 

0744 

1 

6 

B 

3 

2 

0801 

32 

36 

A 

2 

2 

0809 

6 

36 

B 

1 

2 

0873 

20 

34 

B 

1 

2 

0876 

18 

27 

C 

1 

2 

0879 

8 

16 

C 

1 

2 

0826 

1 

?(S) 

0916 

4 

B  S 

2 

2 

0903 

1 

0901 

13 

A 

1 

2 

0909 

15 

51 

B 

1 

2 

0944 

18 

33 

B 

3 

2 

0961 

11 

11 

B 

1 

2 

0964 

12 

16 

D 

1 

2 

0965 

8 

13 

B 

2 

2 

0973 

73 

152 

C 

1 

2 

0976 

34 

49 

C 

1 

2 

0979 

47 

77 

C 

1 

2 

0973 

1 

?(S) 

0972 

1 

?(S) 

1001 

3 

8 

A 

2 

2 

1009 

18 

37 

B 

1 

2 
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LIST  OF  WORKING  PLACES  -  (Continued) 


Working 

Place 

No. 

No.  of 

Jan. / 63 
June/ 64 

Cont.  Rec. 

Feb. / 62 

Nov . / 64 

Mining 

Method 

Ore 

Control 

Ground 

Control 

Breaking 

1101 

28 

51 

A 

2 

2 

1109 

7 

17 

B 

1 

2 

1201 

8 

A 

2  &  3 

2 

1209 

36 

57 

B 

3 

2 

1301 

1 

12 

A 

2 

2 

1309 

30 

42 

B 

3 

2 

1401 

5 

5 

A 

1 

2 

2 

1409 

1 

B 

1 

1 

2 

660 


1,275 
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TABLE  3 

KEY  TO  COMPUTER  OUTPUT  OF  MULTIPLE  REGRESSION  PROGRAMS 


(Program  60  Appendix  2  used  as  an  Example) 


Page  1 

Means,  Variances  and  Standard  Deviations 

Line  1  -  Col  1  -  Program  No.  =  60 

Col  2  -  Variable  No.  =  1 

Col  3  -  Mean  Variable  No.  1 

Col  4  -  Variance  of  Variable  No.  1 

Col  5  -  Standard  deviation  of 
Variable  No.  1 


0.88730290  E  01 

8.8730290 

779.32906 

27.916466 


Pages  2-4 

Simple  Correlation  Coefficients 
Page  2 

Line  1  -  Col  1  -  Program  No.  =60 

_  Variable  No.  notation 

Col  3) 


Col  4  -  ri2) 
Col  5  -  r^) 
Col  6  -  r-j^) 
Col  7  -  ri5) 
Co!  8  -  r16) 
Col  9  -  r^y) 


Where  ^ .  indicates  the  simple 
correlation  of  variable  "i"  with 
variable  "j" 


Pages  5-24  Regression  Coefficients  and  Analysis  of  Variance 
Page  5  -  The  bottom  11  lines  -  Regression  Step  4 


Col  1 
Col  2 

Col  3 


Program  No.  =  60 

Step  No.  =  4  -  Regression  Statistics  - 

Step  4 

Lines  2,  3,  4  and  5  -  Order  of  entry  of 

variables  into  the 
regression  =  1,2, 3, 4. 
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TABLE  3  (continued) 


(program  60  Appendix  2  used  as  an  Example 


Page  5  -  continued 


Col  4  -  lines  2,  3S  4,  and  5  -  No.  of  the  variables 
which  were  entered  into  the  regression  1st, 
2nd,  3rd,  and  4th 


Col  5 


line  1  -  Estimated  coefficient 


line  2  -  Estimated  coefficient 


line  3  -  Estimated  coefficient 


line  4  -  Estimated  coefficient 


line  5  -  Estimated  coefficient 


=  $795.96090 
of  Xl3,  b13 

=  1.0389676 

of  Xl5  bx 

=  6.7959579 

of  Xg,  bg 

-  0.75326410 

of  X20,  b2Q 

=  1.1705417 


Col  6  -  lines  2  -  5  -  Standard  errors  of  estimate 
of  the  coefficients  in  Col  5. 

Col  7  -  lines  2  to  5  =  Col  5  -  Col  6  =  t  ... 

a  -  13,  1,  8,  20  a(4) 

Col  8  -  lines  2  to  5  =  Percentage  variance  reduction 
due  to  the  addition  of  variable  13,  1,  8  &  20 
reduction  due  to  the  addition  of  variables 
13,  1,  8  and  20 

-  line  6  =  Total  of  the  percentage  variance 
reductions . 
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TABLE  3  (continued) 


(Program  60  Appendix  2  used  as  an  Example 


Page  5  -  continued 


Analysis  of  Variance  -  Step  4 


Col  4  -  line  7  =  No.  of  parameters  in  regression 

=  P  =  4 

=  Degrees  of  freedom  for  regression 

-  line  8  =  Degrees  of  freedom  for  error 

=  N-P-l 
=  1270 

(N  =  number  of  observations) 

-  line  9  =  Total  degrees  of  freedom 

=  N-l 
-  1274 


Col  5  -  line  7  =  Sum  of  squares  for  regression 

=  365,779,400.0 

-  line  8  -  Sums  of  squares  for  error 

=  215,208,180.0 

-  line  9  -  Total  sums  of  squares 

=  580,987,580.0 


-  line  10  -  The  standard  error  of  the  dependent 

variable 

=  sn  =  $411.64944 

-  line  11-  The  determinant  of  the  submatrix 

formed  by  simple  correlation 
coefficients  of  the  variables 
entered  in  the  regression 
(13,  1,  8  and  20) 


-  The  estimated  mean  square  for 
regression 

=  91,444,850.0  (Col  5  *  Col  4) 


col  6  -  line  7 
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TABLE  3  (continued) 


(Program  60  Appendix  2  used  as  an  Example 
Page  5  -  continued 

Col  6  -  line  8  -  The  estimated  mean  square  for  error 

=  169,455.26  (Col  5  4  Col  4) 

-  line  10  -  The  sum  of  squares  for  regression 

in  step  4  less  the  sum  of  squares 
for  regression  in  the  preceding  step 
-  20,435,684.0 

=  365,779,400.  -  345,343,720. 

Col  7  -  line  7  -  "F4"  ratio  =  539.64 

-  line  10  -  8,t4"  -  120.60 

Col  8  -  line  7  -  Probability  that  "F1"  could  be  equal 

to  zero=0. 00000 

Page  25  -  Summary  Statistics  on  Residual 

All  statistics  are  identified  on  the  left  except  the 

“'Correlation  of  Residuals'1  which  are  in  the  order 

rXy,  rxz,  ryz,  and  rxz  where  x,  y  and  z  refer  to  the  sequences 

x  -  (Ri>  r2s  *  *  •  *  »  Rn~2^ 

y  —  (R-2 9  R2  9  •  •  •  •  >  Rn~l^ 

z  —  (R3 >  &4 »  •  •  •  •  s  &n) 


APPENDIX  4 


COVARIANCE  ANALYSIS  OF  MULTIPLE  REGRESSION  ON 
GLASSIFICATIONS  OF  STOPE  CONTRACT  DATA 
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APPENDIX  4 

COVARIANCE  ANALYSIS  OF  MULTIPLE  REGRESSION  ON 

CLASSIFICATIONS  OF  STQPE  CONTRACT  DATA 


(Reference  -  Bennett  &  Franklin,  "Statistical  Analysis 
in  Chemistry  and  the  Chemical  Industry"  pp.  457  -  461) 


GIVEN :  Stope  contract  data  representing  a  number  of  fixed 

measurements,  X.  .  on  a  variate  Y,  have  been  classi- 

ij  i 

fied  into  "p"  classifications.  Multiple  regressions 
have  been  performed  on  the  individual  classifications 
and  the  full  set.  The  residual  sum  of  squares  at  the 
completion  of  the  full  regression  has  been  calculated 
in  all  cases. 


REQUIRED: 


METHOD : 


Find  if,  on  the  whole,  the  classifications  differ 

significantly  from  the  full  set  as  regards  to  the 

variability  remaining  after  the  same  parameter 

multiple  regression  has  been  performed. 

Analysis  of  Covariance 
Multiple  Regression 

Source  of  Estimate  Corrected  Sum  of  Squares  D.F. 

P  -  1 

N  -  n  -  p 


Between  Classes 


Within  Classes 


Rn 


&  r- 

P  P 


R 


Total 
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METHOD  -  (continued 

Where  R 

P 


N 


residual  sum  of  squares  after  fitting  of  "n" 

th 

parameter  regression  function  to  the  p 
classification. 

residual  sum  of  squares  after  fitting  of  "n" 

classification. 

total  number  of  observations 


TABLE  6 

ANALYSIS  OF  COVARIANCE 

TIME  SERIES  CLASSIFICATION  (PROGRAMS  81  -  87,  60) 


Source  of  Variation  Corrected  Sums  of  Squares  D  F  MS 

Between  Classes  14,789,877.3  6  2,464,969.6 

i 

Within  Classes  75,560,630.7  1236  61,133.2 

Total  90,350,448.0 


F  (6,1236)  =  40.32 
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TABLE  7 

ANALYSIS  OF  COVARIANCE 
MINING  METHOD  CLASSIFICATION 
(Programs  60  -  67) 


Source  of  Variation  Corrected  Sums  of  Squares  D  F  MS 

Between  Classes  20,588,849.4  6  3,431,474.8 

Within  Classes  69,761,598.6  1236  56,441.4 

Total  90,350,448.0 


F  (6,1236)  =  60.80 


TABLE  8 

ANALYSIS  OF  COVARIANCE 
INCENTIVE  RATE  CLASSIFICATION 
(Programs  60,  70  -  77) 


Source  of  Variation 

Corrected  Sums  of  Squares 

D  F 

MS 

Between  Classes 

75,179,317.7 

7 

10,739,902 

Within  Glasses 

15,171,130.3 

1235 

12,284 

Total 

90,350,448.0 

F  (7,1225)  =  8,743.0 


